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O I We summarize the lessons learned from studies of hard scattering processes in high-energy 

^ ' electron-proton collisions at HERA and antiproton-proton collisions at the Tevatron, with 

the aim of predicting new strong interaction phenomena observable in next— generation 
experiments at the Large Hadron Collider (LHC). Processes reviewed include inclusive 
deep-inelastic scattering (DIS) at small x, exclusive and diffract! ve processes in DIS and 
hadron-hadron scattering, as well as color transparency and nuclear shadowing effects. A 
unified treatment of these processes is outlined, based on factorization theorems of quan- 
^ ' tum chromodynamics, and using the correspondence between the "parton" picture in the 

■ infinite-momentum frame and the "dipole" picture of high-energy processes in the target 

rest frame. The crucial role of the three-dimensional quark and gluon structure of the nu- 
cleon is emphasized. A new dynamical effect predicted at high energies is the unitarity, or 
black disk, limit (BDL) in the interaction of small dipoles with hadronic matter, due to the 
increase of the gluon density at small x. This effect is marginally visible in diffract! ve DIS 
at HERA and will lead to the complete disappearance of Bjorken scaling at higher energies. 
In hadron-hadron scattering at LHC energies and beyond (cosmic ray physics), the BDL 
will be a standard feature of the dynamics, with implications for (a) hadron production 
at forward and central rapidities in central proton-proton and proton-nucleus collisions, 
in particular events with heavy particle production (Higgs), (b) proton-proton elastic scat- 
tering, (c) heavy-ion collisions. We also outline the possibilities for studies of diffractive 
processes and photon-induced reactions (ultraperipheral collisions) at LHC, as well as pos- 
sible measurements with a future electron-ion collider. 
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1 Introduction 



In understanding the nature of strong interactions, progress has mostly come from the investiga- 
tion of certain "extreme" kinematic regions, in which the dynamics simplifies. One such region 
is high-energy hadron-hadron scattering, in which the center-of-mass energy is significantly 
larger than the masses of the hadronic systems in the initial and final state. Historically, this 
was the first area in which powerful mathematical methods, such as dispersion relations and 
Reggeon calculus, could be applied to strong interaction phenomena. They are based on the 
general principles of unitarity of the scattering matrix (conservation of probability) and ana- 
lyticity of scattering amplitudes (causality). These methods have given us important insights 
into general properties of high-energy processes, such as the increase of the radius of inter- 
action with energy predicted by V. Gribov ID El, the Froissart bound for the growth of total 
hadronic cross sections with energy Q, and the Pomeranchuk theorem of asymptotic equality 
of particle-particle and particle-antiparticle cross sections. 

Further progress came with the study of "hard" scattering processes, characterized by a mo- 
mentum transfer significantly larger than the typical mass scale associated with hadron struc- 
ture, fi (a reasonable numerical value for this scale is the p meson mass). Such processes can 
be described in quantum chromodynamics (QCD), the field theory of interacting quarks and 
gluons, the fundamental property of which is the smallness of the effective coupling constant 
in small space-time intervals (asymptotic freedom) |S1I3- Hard processes happen so "rapidly" 
that they do not significantly change the environment of the interacting quarks and gluons inside 
the hadrons. This allows one to calculate their amplitudes using a technique called factoriza- 
tion — a systematic separation into a hard quark-gluon scattering process and certain functions 
describing the distribution of quarks and gluons in the participating hadrons. The simplest 
such process is deep-inelastic lepton-hadron scattering (DIS) in the so-called Bjorken limit, 

~ W"^ ^ /i^, see Fig.[ll Historically, the observation of scaling behavior in the structure 
functions of inclusive DIS |6| gave the first indication of the presence of quasi-free, pointlike 
constituents in the proton |7 1. Another class of processes for which factorization is possible are 
certain hard processes in hadron-hadron scattering, such as the production of jets with large 
transverse momenta or large-mass dilepton pairs. 

A particularly interesting region of strong interactions are hard scattering processes in the 
region where the center-of-mass energy becomes large compared to the momentum transfer, 
W"^ ^ S> /U^. In DIS this limit corresponds to values of the Bjorken variable x ^ 1 
(see Fig. [l]), whence this field is known as "small-x physics." On one hand, because of the 
large momentum transfer, such processes probe the quark and gluon degrees of freedom of 
QCD. On the other hand, they share many characteristics with high-energy hadron-hadron 
scattering, such as a large spatial extension of the interaction region along the collision axis 
(this will be explained in detail in Sec. 12.21) . The treatment of such processes generally requires 
a combination of the methods of QCD factorization and "pre-QCD" methods of high-energy 
hadron-hadron scattering for modeling the dynamics of the hadronic environment of the quarks 
and gluons participating in the hard process. From the point of view of QCD, the high-energy 
(small-a;) region corresponds to a greatly increased phase space for gluon radiation as compared 
to a; ~ 1. QCD predicts a fast rise of the gluon density in the nucleon with decreasing x, and thus 
a strong increase of the DIS cross section with energy lEHHl. A challenging question, which is 
presently being addressed in different approaches, is the role of unitarity of the scattering matrix 
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Figure 1: The kinematics of deep-inelastic lepton-hadron scattering (DIS). The interaction 
proceeds by exchange of a virtual photon, whose four-momentum is given by the difference 
of the lepton momenta, q = k' — k. The hadronic scattering process is characterized by two 
kinematic invariants, the photon virtuality, = —Q'^ < 0, and the photon-proton center-of- 
mass energy, W, or, alternatively, the Bjorken scaling variable, x. (We neglect the target mass.) 

in such processes at high energies. More generally, one hopes to eventually understand the 
quark-gluon dynamics underlying such general high-energy phenomena as the growth of the 
radius of interaction with energy, and the Froissart bound. This dynamics would correspond to a 
strongly interacting quark-gluon system at small coupling constant, and represent a fascinating 
new form of "QCD matter" which could be produced in the laboratory. 

An intuitive understanding of the dynamics of hard scattering processes and the basis for 
QCD factorization can be developed by following the space-time evolution of the reactions in 
certain reference frames. At high energies (small x), the space-time evolution of DIS can be 
discussed in two complementary ways. In a frame where the proton is fast-moving one obtains 
the well-known parton picture of hard processes, in which the hard scattering process involves 
quarks and gluons carrying a certain fraction (here, x) of the proton's momentum. In the proton 
rest frame, on the other hand, the DIS process takes the form of the scattering of a quark- 
antiquark dipole from the target, with the dipole formed a long time before reaching the target, 
and having a distribution of transverse sizes extending down to values ~ l/Q- This represen- 
tation reveals a close relation between DIS at small x and the so-called "color transparency" 
phenomenon — the transparency of hadronic matter to the propagation of spatially small color- 
singlet configurations, as observed e.g. in the suppression of the interaction of heavy quarkonia 
with hadronic matter. The correspondence between the "parton" and the "dipole" picture of 
small-x processes is a powerful tool for analyzing the dynamics of strong interactions in this 
regime. (A pedagogical introduction to these concepts will be given in Sec.|21) 

The experimental investigation of small-x processes became possible with the advent of 
high-energy colliders (colliding beam facilities). Extensive studies of DIS at small x have 
been performed at the HERA electron-proton (ep) collider at DESY. Measurements of inclusive 
cross sections have spectacularly confirmed the rise of the gluon density in the proton at small 
X, as predicted by QCD, down to values x ~ 10~^. Measurements of exclusive processes in 
DIS, such as heavy and light vector meson production (J/V^, p), provide information about the 
spatial distribution of partons in the transverse plane ("generalized parton distributions") and 
allow us to construct a three-dimensional image of the quark and gluon structure of the nucleon. 
Finally, measurements of diffractive processes in DIS, in which the produced hadronic system 
is separated from the target remnants by a large rapidity gap, allow one to probe the interaction 
of various small-size color-singlet configurations with the proton in much more detail than 
inclusive DIS. 
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Another — potentially much more powerful — laboratory for studying small-x physics are 
high-energy proton-proton (pp) and antipro ton-pro ton (pp) colliders, such as LHC at CERN 
(under construction) and the Tevatron at Fermilab. QCD factorization can be applied to pp/pp 
collisions with hard processes, such as the production of dijets with large transverse momenta or 
heavy particles (W^ bosons, Higgs bosons, etc.), which originate from binary collisions of par- 
tons in the two colliding hadrons. At LHC, such processes can probe parton distributions down 
to values of x ~ 10"''. Even higher energies are reached in collisions of cosmic-ray particles 
near the Greisen-Zatsepin-Kuzmin cutoff |10| with atmospheric nuclei. In pp scattering, as 
compared to ep, one is dealing with collisions of two objects with a complex internal structure. 
This results e.g. in a high probability of multiple hard scattering processes at high energies, and 
a much richer spectrum of soft hadronic interactions. Thus, while QCD factorization can still be 
applied to certain hard processes in pp/pp collisions, the modeling of the hadronic environment 
of the quarks and gluons participating in the hard process (or processes) is generally much more 
challenging than in ep scattering. 

This review is an attempt to summarize what has been learned about small-x physics from 
experiments at HERA and the Tevatron and related theoretical studies, and use this information 
to make predictions for new QCD phenomena observable at LHC. With LHC about to be com- 
missioned, and the HERA program nearing completion, this is a timely exercise. It is not our 
aim to give a comprehensive overview of the existing collider experiments and their numerous 
implications for our understanding of QCD. Rather, we identify certain specific "lessons" which 
are of particular importance in making the transition from HERA to LHC: 

• Black-disk limit in dipole-hadron interactions. Because of the fast rise of the gluon den- 
sity at small x, the strength of interaction of small color-singlet configurations (dipoles) 
with hadronic matter can approach the maximum value allowed by s-channel unitarity. 
We quantify this effect by formulating an optical model of dipole-hadron scattering, in 
which the unitarity limit corresponds to the scattering from a "black disk", whose ra- 
dius increases with energy (black disk limit, or BDL). We argue that the onset of the BDL 
regime can be seen in the diffractive DIS data at the upper end of the HERA energy range, 
as well as in elastic pp/pp scattering at Tevatron energies. In DIS at higher energies, the 
BDL leads to a breakdown of Bjorken scaling. In hadron-hadron collisions at LHC ener- 
gies and beyond (cosmic ray physics), the dynamics will be deep inside the BDL regime, 
with numerous consequences for the hadronic final states. 

• Small transverse area of leading partons. Studies of hard exclusive processes at HERA 
show that partons in the nucleon with x > 10"^ and significant transverse momenta 
are concentrated in a small transverse area, <^ 1 fm^, substantially smaller than the area 
associated with the nucleon in soft (hadronic) interactions at high energies. The resulting 
"two-scale picture" of the transverse structure of the nucleon is essential for modeling 
the hadronic environment of the colliding partons in high-energy pp collisions with hard 
processes. 

Based on these observations, we make several predictions for new strong interaction phenomena 
observable in pp, pA (proton-nucleus), and AA (nucleus-nucleus) collisions at LHC: 
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• Hard processes as a trigger for central pp collisions. In pp scattering at LHC, hard QCD 

processes involving binary collisions of partons with momentum fractions Xi 2 > 10~^ 
occur practically only in pp events with small impact parameters (central collisions). This 
makes it possible to trigger on central pp events by requiring the presence of a hard dijet 
(or double dijet) at small rapidities. 

• Black-disk limit in central pp/pA/AA collisions at LHC. The approach to the BDL at 
high energies will strongly affect the dynamics of central pp/pA/AA collisions at LHC 
energies and above. A crucial point is that in hadron-hadron collisions at such energies 
one is dealing mostly with gluon-gluon dipoles, whose cross section for scattering from 
hadronic matter is 9/4 times larger than that of the quark-antiquark dipoles dominating ep 
scattering. We argue that as a consequence of the approach to the BDL the leading partons 
in central pp collisions will acquire large transverse momenta (p5_ ~ several 10 GeV^) 
and fragment independently, resulting in the disappearance of leading hadrons with small 
transverse momenta at forward/backward rapidities, increased energy loss, and increased 
soft particle production at central rapidities. These observable effects allow for experi- 
mental studies of this fascinating new regime of "strong gluon fields" at LHC. We also 
outline the role of the BDL in heavy-ion collisions and cosmic ray physics. 

• Dijfraction in high-energy pp collisions. LHC offers the possibility to study a wide vari- 
ety of diffractive processes in high-energy pp scattering, which probe the interaction of 
small-size color singlets with hadronic matter and can be used to map the gluon distri- 
bution in the proton. Such processes involve a delicate interplay between hard (partonic) 
and soft (hadronic) interactions. A crucial ingredient in understanding the dynamics is 
the information about the transverse spatial distribution of gluons obtained from exclu- 
sive vector meson production in DIS at HERA. 

We also comment on the potential of LHC for parton distribution measurements at small x, and 
for studies of small-x dynamics via photon-induced reactions in ultraperipheral pp/pA/AA 
collisions. Finally, we discuss the opportunities for studies of small-x dynamics provided by 
the planned electron-ion collider. 

The primary purpose of LHC is the search for new heavy particles (Higgs bosons, super- 
symmetry) in high-energy pp collisions. The small-x phenomena we describe here directly 
impact on this program. Heavy particles are produced in hard partonic collisions. For the 
reason described above, the production of heavy particles in inelastic pp collisions happens pre- 
dominantly in central collisions, which are strongly affected by the approach to the BDL, and 
the strong interaction background may be completely different from what one would expect 
based on the naive extrapolation of existing data (Tevatron). Likewise, the search for Higgs 
bosons in diffractive pp events depends crucially on the understanding of the strong interaction 
dynamics in these processes. 
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2 QCD factorization and the space-time evolution of small-x 
scattering 



2.1 QCD factorization of hard processes 

We begin by introducing the basic concept of QCD factorization of hard processes, and outlining 
the space-time evolution of small-s scattering processes at high energies. Our main point is the 
correspondence between the "parton" picture of hard processes in a frame in which the nucleon 
is moving fast, and the "dipole" picture of high-energy processes in the target rest frame. In 
this Section we illustrate this correspondence using as an example the simplest high-energy 
process, inclusive DIS at small x. Below we shall apply these results to analyze the HERA DIS 
data (Sec. 13), and generalize them to processes with exclusive (Sec. HI) and diffractive (Sec. 13) 
final states. The correspondence between the two pictures plays a crucial role in formulating 
the the approach the unitariy limit at high energies (Sec.|6l), and for understanding the dynamics 
of high-energy hadron-hadron collisions (Sees. IT] and El). 

DIS is essentially the scattering of a virtual photon (7*) from a hadronic target, see Fig.[ll By 
the optical theorem of quantum mechanics, the total 7*^ cross section is given by the imaginary 
part of the forward scattering amplitude (virtual Compton amplitude), see Fig.EK- We consider 
the Bjorken limit, in which both the photon virtuality and the 7*]? center-of-mass energy be- 
come large compared to the typical hadronic mass scale, ~ W"^ ^ jj?. As a consequence 
of the asymptotic freedom of QCD, DIS in this limit can be described as the scattering of the 
virtual photon from quasi-free quarks (and antiquarks) in the proton. In the simplest approxi- 
mation, one neglects the interactions of the quarks altogether, and considers the scattering from 
a free quark, see Fig. Eh. This is equivalent to the space-time picture of DIS expressed in the 
parton model [T| . Its basic assumption is that, in a reference frame where the proton moves with 
a large velocity, the interaction of the 7* with the quarks ("partons") is instantaneous compared 
to the characteristic time of their internal motion in the proton. In this picture, the total cross 
section for 7*]? scattering in the Bjorken limit is given by 

A 2 

^rP^X(g2^^) = g^^^^2(x), F,{x) = ej X [g^ (x) + (x)] , (1) 

where a^m is the electromagnetic fine structure constant. Here, e/ are the quark charges (/ = 
u,d,s... labels the quark flavor), and g/(x) and g/(x) are the parton densities, describing the 
number density of quarks and antiquarks carrying a fraction, x, of the fast-moving proton's 
momentum. The transverse momenta of the quarks and antiquarks are of the order kj_ ~ yU^, 
and are integrated over. Eq. ([T]) exhibits the famous property of Bjorken scaling [6], i.e., the 
structure function, F2, depends on the kinematic invariants characterizing the initial state only 
through the dimensionless Bjorken variable (we neglect the nucleon mass). 

This is a direct consequence of the scattering from pointlike, quasi-free particles. 
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(c) (d) 

Figure 2: (a) The total cross section for 7*p scattering is given by the imaginary part of the 
forwad scattering amplitude, (b) The total cross section in the parton model (Bjorken scaling). 
(c,d) QCD radiative corrections, giving rise to the leading scaling violations in as \n{Q'^/Ql). 

The parton model assumption of widely different timescales for the 7*-quark interaction 
and the internal motion of the quarks in the proton becomes invalid in quantum field theory, 
where the ultraviolet divergences introduce a scale larger than Q^. At its most elementary, this 
is the reason why Bjorken scaling is violated in quantum chromodynamics — an argument 
originally due to V. Gribov. The leading scaling violations in ln(Q^/(5o) arise from gluon 
bremsstrahlung, as described by the ladder-type Feynman diagrams shown in Fig. |2j; and d, 
and can be summed up in closed form. Here, as is the strong coupling constant, and is 
an arbitrarily chosen initial scale in the region of approximate Bjorken scaling. The result can 
be expressed in the form of a Q^-dependence of the parton densities, governed by a differ- 
ential equation, the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation 
ifTTl fT2l l9ll. This formulation allows one to retain the basic space-time picture of the parton 
model while incorporating QCD radiative corrections by way of a Q^-dependence of the parton 
densities. Note that through evolution the gluon density in the proton effectively enters into the 
structure functions of 7*p scattering, see Fig. Ell. 

The basic structure of the 7*p total cross section in the Bjorken limit in QCD is that of 
a product of a "hard" photon-parton cross section (involving virtualities ~ Q^) and a "soft" 
matrix element (involving virtualities ~ fj?), describing the distribution of partons in the pro- 
ton. QCD radiative corrections can be incorporated by a systematic redefinition of the "hard" 
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and "soft" factors. This property is referred to as factorization. Factorization for the 7*p total 
cross section in QCD has been formally demonstrated in several different approaches, includ- 
ing operator methods in which the parton densities appear as nucleon matrix elements of certain 
non-local light-ray operators, and the evolution equations coincide with the QCD renormaliza- 
tion group equations for these operators. The calculations have also been extended to sum up 
next-to-leading (NLO) corrections in a., ln(Q^/(5o)- We shall see below that the basic tech- 
nique of factorization can be applied also to exclusive (Sec.|H) and diffractive (Sec.|51) final states 
in 7*p scattering, as well as to certain hard processes in hadron-hadron scattering (Sec.|71). 

More generally, QCD factorization allows one to perform an asymptotic expansion of the 
DIS structure functions in the Bjorken limit. The contribution from the diagram of Fig. |2b, 
Eq. (HJ), determines the leading power behavior at large Q^, with an additional logarithmic 
dependence appearing due to radiative corrections. Fig |2t and d. In the context of operator 
methods this is known as the leading-twist approximation. Power corrections of the order 
^ IQ^ (higher-twist corrections) arise from taking into account the effect of the quark trans- 
verse momentum on the hard scattering process and the interaction of the "struck" quark with 
the non-perturbative gluon field in the proton; the two effects are intimately related because of 
gauge invariance in QCD [||T3|. 

From a mathematical perspective, Bjorken scaling of the moments of the DIS structure func- 
tions can be seen as a consequence of the conformal invariance of the QCD Lagrangian. The 
ultraviolet divergences associated with radiative corrections give rise to anomalous representa- 
tions of the conformal group, with a logarithmic scale dependence. Later we shall see that at 
high energies a new dynamical scale appears in QCD, related to the gluon density in the nucleon 
and its transverse area, which breaks the conformal invariance, and thus leads to the complete 
disappearance of Bjorken scaling — the black disk limit (BDL), see Sec. 

2.2 Space-time evolution of small-x scattering in the target rest frame 

We now turn to DIS at high energies, W'^ ^ ^ ji^, which corresponds to values of the 
Bjorken variable x ^ 1. While this processes can be discussed within the standard QCD 
factorization approach described above, one faces the practical question at which point higher- 
twist (1/Q^-) corrections enhanced at small x, or radiative corrections beyond the DGLAP 
approximation giving rise to factors ln(l/x), become important. These and other questions 
can be addressed in a transparent way by considering the time evolution of DIS in the target 
rest frame, where the process takes the form of the scattering of a small-size qq dipole from 
a hadronic target. More generally, this formulation suggests a new understanding of QCD 
factorization, closely related to the so-called "color transparency" phenomenon observed in 
diffractive processes in hadron-hadron scattering. 

QCD factorization in DIS and the DGLAP approximation have been formulated using the 
covariant language of Feynman diagrams. A typical Feynman diagram relevant at small x is 
shown in Fig. In order to arrive at a space-time interpretation one needs to perform the 
integration over the "energy" variable using the residue theorem. It is this step which actually 
introduces the dependence of the amplitudes on the reference frame. Alternatively, one may 
directly trace the space-time evolution using the language of time-ordered perturbation theory. 
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Figure 3: (a) A typical Feynman diagram for inclusive scattering at small x. (b) Time 
evolution of small-x scattering in the target rest frame, (c) The decay of a fast parton in the 
qq dipole. The degradation of longitudinal momenta is accompanied by a random walk in the 
transverse coordinate. 

In this formulation, the time scales of the processes are determined by the energy denominators 
associated with the various transitions, via the energy-time uncertainty relation. At = 1/A£J. 

In the target rest frame, the virtual photon in DIS at small x moves with 3-momentum 
P Q"^ /(2mNx). Consider its conversion into a qq pair with longitudinal momenta zP and 
(1 — z)P and transverse momenta ±k± ("longitudinal" and "transverse" are defined relative to 
the direction of motion of the photon). The energy denominator for this transition is 

= (3) 

where M^^ = (m^ + kjj/lz^l — z)] is the invariant mass of the qq pair (m^ is the quark mass, 
k± = \k±\). For a longitudinally polarized photon, the dominant contribution to the cross 
section comes from values z ~ 1/2 and k± ~ Q, for which M|g ~ (this will be explained 
in more detail in Sec. 13.21) . For such values the time associated with the 7* qq transition is 

At = 1/AE ~ Q^/P ^ l/{2mNx). (4) 

At small X, the photon converts into a qq pair long before reaching the target, as illustrated in 
Fig. 13). Both the quark and antiquark move essentially with the speed of light. The distance 
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between the point of their creation and the target, the so-called coherence length, is given by 



Icob = cAt. (5) 

It is important to realize that the qq wave packet remains well localized in the longitudinal 
direction as it travels towards the target, and that the transverse separation between the quark 
and antiquark (at a given time) is a meaningful concept. It is generally of the order d ~ ^/Mgg ~ 
1/Q, and thus small if is sufficiently large. 

In short, high-energy 72-hadron scattering in the target rest frame is essentially the scatter- 
ing of a small-size qq dipole from the target hadron. In this formulation, the quantity describing 
the strong interaction effects is the cross section for dipole-hadron scattering. It can be com- 
puted using methods of QCD factorization, with the dipole size, <^ /i acting as the factor- 
ization scale, see Sec. 13.21 Without going into details, we can immediately state an important 
property of the cross section, namely that in the limit d it vanishes as 

up to logarithmic corrections in d. Eq. Q reflects a fundamental property of QCD as a gauge 
theory — the interaction of a small-size color singlet object with hadronic matter is small 
("color transparency"). In this understanding of QCD factorization, high-energy 7*-hadron 
scattering exhibits a close relation to the interaction of heavy quarkonia with hadronic matter 
and a number of other color transparency phenomena in hadron-hadron scattering. 

To make the dipole picture quantitative, one has to take into account the effects of QCD 
radiation. In particular, this is necessary in order to determine the coefficient in Eq. ^ with 
logarithmic accuracy. The importance of different types of radiation can again be studied using 
the language of time-ordered perturbation theory in the target rest frame. The characteristic 
time for the quark to radiate a gluon with longitudinal momentum fraction Xg and transverse 
momentum k±^g, relative to the time the qq pair spends between its creation and "hitting" the 
target, ©, is (cf. Fig. Ob) 



Ati _ r„o m^ + fc^H 



2 I g 



zil 



Xg Xg 



(7) 



If X is sufficiently small, and for average values of z, the emission process can be repeated 
several times before the evolved system reaches the target. 

There exist several kinematic domains where gluon emission during the propagation of the 
qq wave packet is likely because of a large phase space at small x. One is the emission of partons 
with transverse momenta smaller than k^ of the parent parton. Each such emission contributes 
a factor as ln(Q^/Qo) '^he amplitude, where the logarithm arises from the integration over 
the phase volume of the radiated gluon. In the standard QCD description of DIS, these are the 
radiative corrections summed up by the DGLAP evolution equations for the parton distributions 
described above ||Tl][T2l|9l (see also Ref. lO). In the context of the dipole picture at small x, 
the summation of these corrections in leading order (LO) corresponds to a dipole-hadron cross 
section of the form lfT5lfT6nl7J 

^a,--hadro„^^^^2) ^ ^FU'as{Qln)xG{x,Ql^). (8) 
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Here = 4/3 is the Casimir operator of the fundamental representation of the SU{3) gauge 
group. Furthermore, as{Ql^) is the LO running coupling constant and G{x, Qeff) the LO gluon 
density in the target. They are evaluated at a scale Ql^ oc d^'^. The coefficient of proportionality 
is not fixed within the LO approximation, and needs to be determined from NLO calculations 
or from phenomenological considerations, see Sections I3^ andl4l' 

Eq. ^ actually quotes a simplified expression for the dipole-hadron cross section. The 
original expression involves an integral over the gluon momentum fractions, which is concen- 
trated in a narrow range above x. Also neglected in Eq. dH) is the contribution proportional to 
the quark/antiquark distribution in the target, which at small x is suppressed compared to the 
gluon distribution. This contribution would lead to a flavor dependence of the dipole-nucleon 
cross section [19 J. 

Another large contribution, specific to small x, comes from the large phase space in rapidity 
(oc Zg) for emission of gluons without strong degrading of transverse momenta in the leading 
approximation. Such emissions give rise to factors as{Nc/2TT)Ay, where N^. = Sis the number 
of colors in QCD, and Ay = {yi — yi+i) is the difference in rapidities between successive 
partons in the ladder. In terms of x, this corresponds to corrections proportional to ln(xo/a:;), 
where xq ~ 0.1 accounts for the fact that nucleon fragmentation enters in the definition of the 
gluon density in the nucleon and does not produce a logarithm in x. If the rapidity interval for 
emissions {i.e., the lifetime of the quark-gluon system) becomes very large, one needs to sum 
these logarithms in addition to the hi{Q'^ /Ql) terms lEUIl . see Sec. l3.3[ 

QCD radiation generally leads to an increase of the transverse size of the "dressed" dipole 
with decreasing x, and thus to an increase of the radius of the dipole-hadron interaction with 
energy. Each individual emission shifts the transverse coordinate of the radiating parton by 
Ap ~ l/fc_L, see Fig. 13:. If there are n successive emissions with comparable, randomly 
oriented fc^ (this is the case in the limit of large In x), the overall shift is \li 

Ap2 = n/kl = yl{Aykl), (9) 

where y is the rapidity of the initial parton. A similar diffusion mechanism for soft partons was 
discussed by V. Gribov as a model for the increase of the radius of soft hadronic interactions 
with energy EIZI- In the case of hard processes such as 7*-hadron scattering, in the region 
where the DGLAP approximation is valid, the rate of expansion with energy is much smaller 
than for soft interactions, because of the larger transverse momenta of the emitted partons and 
the larger rapidity intervals between the emissions. This manifests itself e.g. in a much weaker 
energy dependence of the t-slope of hard exclusive processes as compared to elastic hadron- 
hadron scattering ||2T1 . see Sec.|4l 

For transversely polarized virtual photons the space-time picture of the interaction is more 
complicated than in the longitudinal case. Owing to the different spin structure of the 7^ qq 
vertex, configurations of very different size — from hadronic size to 1/Q — contribute to the 
interaction. The hadronic size configurations correspond to 2 ~ 1 or 0, and k±_ ~ Aqcd- They 

'There is an approach to high-energy scattering in which the projectile particle is represented as a superposition 
of eigenstates of the scattering matrix, see e.g. Ref. 1 18J and references therein. Equation (|8j implies that states 
with different transverse size, d, should be orthogonal. However, the extension of Eq. [S] to the case of elastic 
scattering indicates that transitions between configurations with different d are allowed for finite t. This suggests 
that the the eigenstate model should be a reasonable approximation only for small values of t. 
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are dual to two jets aligned along the virtual photon direction and are referred to as aligned 
jet configurations. They are expected to interact with the target with typical hadronic cross 
sections, giving the dominant contribution to the structure function F2 at ~ few GeV^ and 
X ~ 10^^, see Sec. 13.21 Also, such configurations can easily scatter elastically from the target, 
and thus are an important source of diffractive scattering, see Sec.|5l 



3 Inclusive Yp scattering at small x 
3.1 DGLAP evolution and the HERA data 

We start our discussion of 7*p scattering at high energies with inclusive DIS. Inclusive DIS is the 
main source of information about the parton distributions in the nucleon at small x. Because 
of the relatively simple structure of QCD factorization, it is also the main testing ground for 
higher-order QCD calculations and resummation approaches. 

The validity of QCD factorization and DGLAP evolution for inclusive DIS have extensively 
been tested in fixed-target experiments, probing the quark/antiquark densities in the nucleon at 
values X > 10^^ see e.g. Ref. ll22ll for a review. Going to smaller x, DGLAP evolution produces 
a fast increase of the parton densities, related to the fact that the gluon has spin 1 JEHU, which 
implies a fast increase of the DIS cross section with energy. This prediction has spectacularly 
been confirmed by the measurements with the HERA ep collider. Fig.|4|shows a summary of the 
F2 proton structure function data taken by HI and ZEUS compared to a QCD fit based on NLO 
DGLAP evolution ll23l . The data clearly support the interrelation of the x- and Q^-dependence 
as predicted by DGLAP evolution. The analysis of the data found that effects of next-to-next- 
to-leading order (NNLO) terms of the form of multiplied by a function of as ln{Q'^/Ql) 
generally appear to be small. It is remarkable that the DGLAP approximation, which does 
not account for all potentially large terms containing ln(l/a;), describes the presently available 
high-energy data so well. 

More detailed insights into the "workings" of the DGLAP approximation can be gained by 
studying the effective power behavior in x of the structure function and the individual parton 
distributions in the NLO fit, 

F2 oc x~^% xG{x) oc x~^\ ^ejxQfix) oc x'^" (x < 10"^), (10) 

/ 

where the exponents depend on Q^, see Fig.|5| At low Q'^, X2 ~ 0.1, reflecting the energy depen- 
dence expected for the cross section of soft hadronic processes. Starting from ^ 0.5 GeV^ 
A2 grows, reaching a value of ~ 0.4 at ~ lOGeV^ (A.Levy, private communication). For 
> 3 GeV^, one observes that A^, ~ A2, indicating that in this Q^-region the x-dependence 
of the structure function is indeed driven by the gluon distribution. For lower Q^, however, Xg is 
significantly different from A2, becoming even negative at ~ 2 GeV^. Thus, while the NLO 
DGLAP approximation formally describes the x-dependence of the structure function even at 
low Q"^, the price to be paid is the lack of a smooth matching of the x-dependence of the gluon 
distribution to the soft regime. This may indicate the presence of significant corrections to the 
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Figure 4: The proton structure function, ^2(0;), as measured by the HI and ZEUS experiments 
at HERA |23 1. Also included are data from fixed-target experiments. The lines show a QCD fit 
based on the NLO DGLAP approximation. 

leading-twist description of DIS at small x for < 3 GeV^. The dynamical origin of these 
corrections will be discussed in Sec. 13.21 

The data show that the deviation from the soft energy dependence of F2 starts at surprisingly 
low scales, <^ 1 GeV^. Within the DGLAP approximation this behavior can be explained 
by the presence of a large non-perturbative gluon density in the nucleon at moderate x at a low 
scale Ii24ll . This is principally consistent with the idea of spontaneous chiral symmetry breaking, 
according to which most of the nucleon mass resides in gluon fields. 
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Figure 5: The exponents characterizing the x-dependence of F2, A2, the gluon distribution, Ag 
(black squares), and the sea quark distributions, Ag (green triangles), cf. Eq. (flUI) . as extracted 
from the NLO DGLAP fit to the HI and ZEUS data (A.Levy, private communication). 

3.2 Space-time picture of inclusive DIS 

Many of the observed features of inclusive DIS at small x can be understood within the space- 
time picture in the target rest frame, cf. Sec. l2.2[ In this formulation, corrections to the leading- 
twist approximation at low Q'^ appear because of the contribution from large dipole sizes. This 
allows us to quantify the region of validity of the leading-twist approximation, and develop an 
"interpolating" approximation valid in a wide range of . 

Following the logic outlined in Sec. 12.21 one can express the total 7*p cross section at small 
X as a superposition of qq dipole cross sections, characterized by the longitudinal momentum 
fraction of the quark, and the dipole size, d: 

(JlA^^Q^) = dz j d^da'''^-^'''^'°''{z,d,x)\tPl^{z,d,Q^)\\ (11) 

where ip1rp(z, d) denotes the light-cone wave function of the qq component of the virtual pho- 
ton, calculable in quantum electrodynamics. An important question is which dipole sizes dom- 
inate in the integral. For a longitudinally polarized photon, the modulus squared of the wave 
function is given by 

|^I(^,rf,g^)r = ^-^^Y.^)[z{l-z)K,{ed)]\ (12) 



15 



where Kq is the modified Bessel function and = z(l — z)Q'^ + mj ll25ll . One can verify by 
direct calculation that in this case the contributions from large dipole sizes are suppressed at 
large Q^, if the integral (fTTTl is evaluated with the LO expression for the dipole-nucleon cross 
section, Eq. (EJ) [21]. In fact, Eqs. (Elinil and (EJ) are formally equivalent to the LO DGLAP 
approximation in QCD, cf. the discussion below. The effective scale in the gluon distribution 
entering the dipole cross section can be determined by comparing (fTTT) with the LO DGLAP 
expression; one finds Qgff ~ 9/(i^ for HERA kinematics [26|. Note that the factor xG{x, Ql^^) 
in Eq. ^ results in a fast increase of the cross section with energy, in contrast to the two-gluon 
exchange model of Refs. I27..28..29J . where the cross section is energy-independent. 

When applying Eq. (fTTT) to transversely polarized photons, the distribution of dipole sizes is 
significantly wider than in the longitudinal case. At Q'^ ~ few GeV^, the transverse cross sec- 
tion receives sizable contributions from dipole sizes for which the perturbative approximation 
for the dipole-nucleon cross section, cf. Eq. ©, becomes invalid. Still, at large the pertur- 
bative contribution should dominate, because of the faster increase with energy of the parton 
distribution for the smaller-size quark-gluon configuration. The contribution from large-size 
qq configurations is strongly suppressed by Sudakov form factors; it is actually represented by 
large-size qqg, . . . configurations. 

Equation (fTTT) can serve as the basis for an "interpolating" model that describes YP inter- 
actions over a wide range of for both transverse and longitudinal polarizations QUI . There 
is ample evidence — e.g. from studies of and vrA^ elastic scattering — that real photons in 
high-energy reactions have transverse sizes comparable to pions. A way to ensure this within 
the qq dipole description is to introduce a dynamical quark mass of ~ 300 MeV, which is 
consistent with the phenomenology of spontaneous chiral symmetry breaking lISTIl. The cross 
section for the scattering of such a "hadronic-size" dipole with the target can then be inferred 
from the vrA^ scattering data. For small dipoles, d < 0.4 fm, the cross section can be calcu- 
lated perturbatively. When evaluating the leading-twist expression, it is important to accurately 
treat the kinematic limits of the integral over the gluon momentum fractions, as this leads to 
an additional dependence of the dipole cross section on Q^. A dipole cross section obtained 
by matching the two prescriptions is shown in Fig.|^ This function is then averaged with the 
photon wave function for massive quarks, cf. Eq. dTTT) . This model reproduces well the HERA 
F2p data for > 0.1 GeV^ and correctly predicts aL ROI . 

In order to make contact with the analysis of Sec. 13.11 we need to state more precisely 
how the dipole picture is related to the DGLAP approximation in QCD. In LO, it has been 
demonstrated explicitly that Eqs. (El) and dTTT) can be obtained by rewriting the LO DGLAP 
expression for the 7*p cross section 1 17|. This simple relation appears because in the leading 
logarithmic approximation the separation of the process according to time in the target rest 
frame — transition of the virtual photon into a qq pair (photon wave function), and interaction 
of the pair with the target (dipole cross section) — coincides with the separation of transverse 
momenta in k^^ ~ Q"^ and A;^ ^ in the partonic ladder. Beyond the leading order, one needs 
to explicitly include qqg component of the photon wave function, and the distinction between 
the wave function and the dipole interaction with the target becomes more delicate. Although 
in principle the leading-twist dipole picture should be equivalent to the DGLAP approximation 
in any order of the expansion, in practice the problem of formulating a consistent dipole picture 
in NLO has not been solved yet. Also, it is worth emphasizing that the correspondence between 
parton distributions and amplitudes of physical processes is not always direct. In particular, the 
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Figure 6: The dipole-nucleon cross section in the "interpolating" model of Ref. f^D"]. Shown are 
the results corresponding to two different values of the parameter, A, determining the effective 
scale, Qeff- 



parton model-type contribution, which naturally leads to diffractive processes, is hidden in the 
boundary condition for the QCD evolution of the parton distributions. 

In spite of the lack of an explicit NLO analysis in the dipole picture, it still seems to be of 
relevance that the region of where the NLO DGLAP analysis leads to a gluon exponent Ag 
dropping below the soft value (see Fig.|5l) corresponds in the dipole model to contributions from 
d > 0.4 fm, where nonperturbative effects become important. Thus, it seems that the "anoma- 
lous" behavior of \g is a consequence of the leading-twist DGLAP approximation trying to 
mock up higher-twist corrections at low Q^. The dipole picture allows us to quantify the region 
of applicability of the leading-twist approximation at low Q^, and suggests a natural way to 
incorporate non-perturbative effects. 

Equations ^ and (fTTT) are valid also within the leading asln(xo/x) approximation. Fur- 
thermore, they can be derived from the eikonal model expression for the propagation of a heavy 
quarkonium through a hadronic medium [32|. 



3.3 Breakdown of the DGLAP approximation at very small x 

The observation of the fast increase of parton densities at small x has stimulated theoretical dis- 
cussions of the stability of the DGLAP approximation at small x. In fact, in the kinematic limit 
of fixed and x the effective parameter of the perturbative QCD expansion is multiplied 
by a factor ln(xo/x), which arises due to gluon emission in multi-Regge kinematics (rapidity 
distance between adjacent gluons ^ 2), and the hierarchy of dominant terms is changed as com- 
pared to the DGLAP approximation. The constant Xq is determined by the typical momentum 
fraction in the initial parton distributions; usually Xq ~ 0.1. 
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A simple kinematic estimate shows that in typical HERA kinematics the DGLAP approx- 
imation is still reliable. The rapidity span at HERA is approximately \i\{Q / xm^) ^ 10 for 
X = 10~^ and Q = 2.5 GeV. To obtain a significant ln(xo/a;) term, the distance in rapidity 
between adjacent partons in the ladder should be ^ 2. Thus, the number of radiated gluons in 
multi-Regge kinematics at HERA is ^ (10 — 4)/2 — 1 ^2, where we took into account that 
each of the fragmentation regions occupies at least two units in rapidity. This simple estimate 
agrees well with a numerical study of NLO QCD evolution, which indicates that the average 
change of x in the HERA region does not exceed 10, corresponding to Ay = \n{xo/x) ~ 2, 
provided that Ql > 1 GeV^ f33\. Since one (two) logarithms of x are effectively taken into ac- 
count by the NLO (NNLO) approximation, there is no need for a special treatment of \n{xo/x) 
effects at HERA kinematics. A similar estimate shows that at LHC kinematics the radiation of 
5-6 gluons is permitted. Thus, at LHC energies and above the resummation of ln(xo/x) terms 
becomes a practical issue. 

The program of resumming leading ctem ln(a;o/a;) terms started in quantum electrodynamics 
lEUI . In QCD, the reggeization of gluons slows down the energy dependence of amplitudes of 
high-energy processes [34, 35|. ^ In the leading ln(xo/a;) approximation (Balitsky-Fadin- 
Kuraev-Lipatov, or BFKL, approximation) [36|, where energy-momentum conservation and 
the running of the coupling constant are neglected, the reggeization of gluons is canceled by 
contributions from multigluon radiation. NLO corrections to a large extent subtract kinemati- 
cally forbidden contributions, leading to a large negative contribution to the structure functions 
ll37l l38l . Another feature of this approach is the lack of an unambiguous separation between 
perturbative and nonperturbative QCD effects |37|. Thus, this approximation seems to be lim- 
ited to the description of single-scale hard processes where DGLAP evolution is unimportant 
in a wide kinematic range, such as YiQ'^) + 7*(Q^) ^ hadrons, or two-body processes where 
the hardness is controlled by proper choice of final state like, such as 7* + 7* — TT. 

The resummation approaches of Refs. ll39l 1401 predict a significantly slower increase of 
amplitudes with energy than the LO BFKL approximation, and possibly even oscillations in the 
energy dependence. Most of the reduction is due to the better account of energy-momentum 
conservation in these approaches, and account of the running of the coupling constant. At 
extremely small x (beyond the reach of LHC) much of the LO BFKL results reappear, but 
with a slower dependence on x. For the parton densities in the nucleon, where a;o ~ 0.1 is a 
reasonable value for the constant in the ln(a;o/a;) factor, resummation effects should be small 
for X > 10~^, that is, for the whole HERA range above > 2 GeV^. At smaller x, the result 
of the resummed evolution is close to that of NLO DGLAP evolution down to x ~ 10^, but 
differs strongly from NNLO [41 1. This suggests that NLO DGLAP evolution could be a good 
guess for the parton densities down to the very small x values probed at LHC, even though the 
underlying dynamics may change significantly at x < 10""^. 

■^The high-energy behavior of two-body ampHtudes with color-octet quantum numbers in the crossed channel 
in QCD is given by the Regge pole formula, {l/x)^'-^\ where P{t) decreases with increase of —t. In leading order 
of as QCD gives P{t) — 1. Thus, gluons in QCD (as well as quarks) are reggeons. 
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4 Exclusive processes in 7*^9 scattering at small x 



4.1 QCD factorization for hard exclusive processes 

The concept of QCD factorization can be extended to certain exclusive channels in 7*p scatter- 
ing, namely processes of the type 

72(g) + N{p) "Meson"(g + A) + "Baryon"(p - A), (13) 

at large virtuality, = — g^, and center-of-mass energy, W'^ = (p + g)^, with fixed x = 
Q'^/(W'^ + Q"^), and fixed small invariant momentum transfer, t = A^. Examples include the 
production of light vector mesons (p, p') [21J, heavy vector mesons (J/V^, ijj', T) ffSll . and real 
photons (deeply-virtual Compton scattering, DVCS) 1021 |43l 031051061071381. Closely related 
to these processes are certain hadron-induced reactions, such as the diffractive dissociation of 
pions, TT + T 2jets + T, where T denotes a hadronic target (nucleon or nucleus) [ 1 6 1 . These 
exclusive processes probe the interaction of small-size color singlets with hadronic matter in 
much more detail than inclusive DIS. They also provide new information about the transverse 
spatial structure of the nucleon, contained in the so-called generalized parton distributions. 

The basis for the analysis of exclusive processes (fT3t is the QCD factorization theorem [49 1, 
which extends the initial analysis of Ref. [21] for the small-x limit. It states that the amplitude 
can be represented as a convolution of three functions, as depicted in Fig.jVl 

^-flN^M+B = / dz dxi fi/p{xi,x - Xi,t;fi) Hij{xi,x,z,Q^;fi) (j)f{z,p) 

+ power corrections. (14) 

Here, / is the generalized parton distribution (GPD), which describes the amplitude for the 
nucleon to "emit" and "absorb" a parton with longitudinal momentum fractions Xi and X2 = 
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Figure 7: Factorization of the amplitude of hard exclusive meson production, Eq. (fT4|l . 
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Xi — X, respectively, accompanied by an invariant momentum transfer, t, and, possibly, a tran- 
sition to another baryonic state. At zero momentum transfer, xi = X2 and t = 0, the GPD's 
coincide with the usual parton densities measured in inclusive DIS. (For a review of the prop- 
erties of GPD's and their applications, see Refs.l SOllSTl .) Furthermore, (p^^ is the distribution 
amplitude describing the conversion of a qq pair with relative longitudinal momentum fraction 
z to the produced meson (or photon). Finally, Hij denotes the amplitude of the hard partonic 
scattering process, which is calculable in powers of ^^(Q^)- The indices i,j label the different 
parton species. The contribution of diagrams in which the hard scattering process involves more 
than the minimum number of partons is suppressed by l/Q^. An important consequence of fac- 
torization is that the t-dependence of the amplitude rests entirely in the GPD. Thus, different 
processes probing the same GPD should exhibit the same t-dependence. 

4.2 Space-time picture of hard exclusive processes 

The physics of hard exclusive processes at small x becomes most transparent when following 
the space-time evolution in the target rest frame. As in the case of inclusive scattering, this 
approach allows one to expose the limits of the leading-twist approximation, and to quantify 
power corrections due to the finite transverse size of the produced meson. 

In exclusive vector meson production, 'jIN VN, one can identify three distinct stages 
in the time evolution in the target rest frame. The virtual photon dissociates into a qq dipole 
of transverse size d ~ l/Q at a time = /coh/c ~ l/(mNx) before interacting with the 
target, cf. Eq. The qq dipole then scatters from the target, and "lives" for a time r/ ^ 
before forming the final state vector meson. The difference in the time scales is due to the 
smaller transverse momenta (virtualities) allowed by the meson wave function as compared to 
the virtual photon. 

In the leading logarithmic approximation in \ia.{Q'^ / Aq^^^^) , the effects of QCD radiation can 
again be absorbed in the amplitude for the scattering of the small-size dipole off the target. It 
can be shown by direct calculation of Feynman diagrams that the leading term for small dipole 
sizes is proportional to the generalized gluon distribution, G{xi,X2, t; Qlff), where Qgff oc d^"^ 
ifTTI . A simpler approach is to infer the result for the imaginary part of the amplitude from 
the expression for the cross section, Eq. via the optical theorem. The imaginary part is 
proportional to the generalized gluon distribution at xi = a; and X2 = 0. At sufficiently large 
Q^, the generalized gluon distribution at small xi and X2 can be calculated by perturbative 
evolution, starting from the "diagonal" generalized gluon distribution, xi = X2 x, at a low 
scale II52I l53l l54l . In applications to vector meson production at HERA, where the effective 
scale is of the order Ql^ ~ few GeV^, the "skewness" effects induced by the evolution are not 
very substantial, and one may approximate the generalized gluon distribution by the diagonal 
one at the scale Ql^^. It is convenient to separate the t-dependence and write the diagonal 
generalized gluon distribution in the form 

G{x,x,t;Qls) = G{x,Qls)Fg{x,t-Q^^s), (15) 

where G(x, Qlf^) is the usual gluon density and Fg is the "two-gluon form factor" of the target, 
which satisfies Fg{x,t = 0; Qeff) = 1- Altogether, one obtains for the dipole-hadron scattering 
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amplitude in this approximation 

A'"^-^{x, d\ t) = 2tiiF^ a.(Qeff) ^^^(x, t, Ql^) Fg{x, t, Ql^). (16) 

The amplitude for the hadronic process ([T3t is then given by the convolution of Eq. ([T6b with 
the light-cone wave function of the virtual photon, Eq. ([T2b . and that of the produced vector 
meson, . In coordinate representation, 

^^*N-.VN ^ dz j d^dijl{z,d) A'''^-^{x,d\t) ij^{z,d), (17) 

where the integration is over the quark longitudinal momentum fraction, z, and the transverse 
dipole size, d. 

In Eq. (flTt . the wave function of the vector meson of transverse size l/my is convoluted 
with the wave function of the virtual photon of significantly smaller transverse size, 1/Q. One 
may say that the meson in this process is "squeezed", i.e., forced to couple in a configuration 
much smaller than its natural hadronic size. In the leading-twist approximation one neglects 
the spatial variation of the vector meson wave function and substitutes it by the distribution 
amplitude, 

^py^z^d) ^ V^^(^,0) = 0^(^) (18) 

(in momentum representation, the distribution amplitude is the integral of the wave function 
over transverse momenta). The integral over transverse sizes can then be performed explicitly, 
using Eqs. (fT2ll and ([T6b . After restoring the real part of the amplitude using its analyticity 
properties, the differential cross section is obtained as [21 J 

dal ^^^^ STT^TvmvVy 

2 

F^{x,t;Ql^). (19) 

Here, Fy is the leptonic width of the vector meson, which defines the normalization of the 
meson wave function, and 

VV ^ I f dz^^^ I \' dzc^^{z), (20) 
2 Jo z(\ -z) I Jo 

?7y — > 1 at asymptotically large Q^. These expressions apply to production by a longitudi- 
nally polarized photon. For transverse polarization, the nonperturbative contribution is sup- 
pressed only by a Sudakov-type form factor, similar to the case of F2{x, Q^) in inclusive 7*p 
scattering. This contribution originates from highly asymmetric qq pairs (z ~ or 1) in the 
7^, wave function, which have transverse size similar to that of hadrons. We note that elastic 
photo/electroproduction of J/ 4' mesons has been evaluated also within the LO BFKL approx- 
imation [55 1. The function xG{x,Q'^) that enters there has no relation to the conventional 
DGLAP gluon distribution, which is defined within the DGLAP approximation only. 

Equation (fT9b is based on the leading logarithmic approximation in ln(Q^/AQ(.j-,), as well 
as on the leading-twist approximation, Eq. dTSt . While it already exhibits many qualitative 
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Figure 8: The average dipole size, d, (left) and the effective scale, Q^^f, (right) in exclusive 
vector meson production (p, J ftp, T) by longitudinally polarized photons, as a function of 
|E^I5^ . Also shown are the shown are the average values of d in the integrand of the expression 
for the inclusive cross section, ctl. 

features seen in the data (see below), two important effects need to be taken into account before 
a quantitative comparison can be attempted. First, because the wave function of the vector 
meson in Eq. (flTt is significantly broader than that of the 72, the effective dipole sizes in 
the meson production amplitude, Eq. (fT6t . are substantially larger than in aL, Eq. (fTTTl . see 
Fig- 111 As a result, the effective scale in the gluon distribution, Q^^f, is smaller in vector meson 
production than in cr^, see Fig. [H |26l |56l|. This effect slows the x- (energy) dependence of 
the cross section compared to the naive estimate, Q^f^ = Q^. Second, numerical studies using 
model wave functions show that retaining the full rf-dependence of the vector meson wave 
function in the convolution integral dTTl) results in a substantial decrease of the absolute cross 
section at moderate as compared to the leading-twist approximation, Eq. (fTSt . as well as 
in a slower dependence ll26ll56l . These higher-twist effects, related to the finite size of the 
vector meson, limit the region of validity of the leading-twist approximation (fTSt and need to 
be taken into account in quantitative estimates at low Q'^. 



4.3 Vector meson production at HERA 

With proper choice of the effective scale, Q^f^, and inclusion of higher-twist effects due to the 
finite transverse size of the meson, one can quantitatively compare the results of the leading 
logarithmic approximation, Eqs. (fT6b and (flTb . with the HERA data on heavy and light vector 
meson production. The data confirm in particular the following predictions of this picture: 

• Increase of cross section with energy. Equation (fT^ implies that da/dt{t = 0) grows 
with energy as [xG{x, Qlff)]^, with Ql^ estimated to be ~ 3 GeV^. When combined with 
the LO gluon density obtained from fits to DIS data, this implies a growth oc Such 
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behavior has been observed for p production at = 10 — 20 GeV^, and for J/^ produc- 
tion starting from = ||57]| . The later onset of the hard regime for p electroproduction 
is due to the rather slow "squeezing" of the qq configuration in the p meson; it reaches 
a size comparable to that of the J/ip only at ~ 20GeV^, see Fig. [H^ The naive 
choice Qeff = would imply a too fast growth, cf. Fig.|5l For soft interactions, on the 
other hand, da /dt(t = 0) oc and the growth is even smaller for the cross section 

integrated over t. 

• Decrease of cross section with Q"^. The decrease with Q'^ of for p-meson production, 
and of the total cross section for J /ip production, is slower than due to the Q^- 
dependence of a^G in Eq. (fT^ . as well as finite-size (higher-twist) effects. This is best 
observed in J /ip electroproduction, where the model of Ref. |55|, which neglects finite- 
size effects, predicts a decrease of the cross section by a factor of ~ 5 faster than observed 
in the kinematic region covered by the HI experiment. 

• Absolute cross sections. The absolute cross sections for vector meson production are well 
reproduced, provided that higher-twist effects due to the finite size of the vector meson 
are taken into account ll26ll56l . 

• Dominance of longitudinal cross section. The data on p production indicate ^ a-p for 

^ rriy, in agreement with our picture. 

• Universality of t-dependence. Comparison of p and J electroproduction data clearly 
show the universality of the t-dependence at large Q^, where the vector mesons are 
"squeezed", and the t-dependence originates solely from the two-gluon form factor, see 
Fig.13 

• Flavor symmetry. Since the interaction of the qq dipole with the gluon distribution is 
flavor blind, one expects the restoration of SU (3) flavor symmetry in vector meson pro- 
duction for Q"^ ^ my. For example, : p = (2 : 9) in the flavor symmetry limit. The 
violation of SU{?>) flavor symmetry due to increase of the wave function of the vector 
meson at small distances with increasing quark mass leads to an enhancement of this 
ratio by a factor ~ 1.2. 

A new situation is encountered in the photoproduction of T mesons. In this case, the 
approximation of the generalized gluon distribution by the usual gluon density becomes in- 
valid (large "skewness" and large Qeff)' ^^id the real part of the amplitude becomes significant. 
Together, these effects increase the predicted cross section by a factor of about 4 ll59l l60ll . 
For T production Qeff ~ 40GeV^, leading to an energy dependence of the cross section as 
da/dt{t = 0) oc W^-\ 

Closely related to vector meson production is the production of real photons (deeply virtual 
Compton scattering, DVCS). This process has been the subject of intense theoretical study in the 
region of moderate x, accessible in fixed-target experiments (HERMES at DESY, COMR^SS 

^In the case of the p— meson production initiated by the transverse photon, the squeezing is generated by the 
Sudakov form factor as well as by the more rapid increase with energy of the small size contribution. The observed 
behavior of cjl /ctt can be fitted within the current models 1581 . 
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Figure 9: The HERA HI and ZEUS data for the t-slopes of the differential cross sections for 
the exclusive electroproduction of p, (p and J/ijj mesons, as a function of Q^. The convergence 
of the different slopes at large Q"^ indicates the dominance of small-size configurations in the 
production process ("squeezing"). The solid line shows the Q ^-dependence obtained in the 
calculation of Ref. ll26l . The data are from Ref. iBTll . 

at CERN, Jefferson Lab), and is considered the main tool for probing the generalized quark dis- 
tributions in the nucleon At small x, the DVCS amplitude has been computed 
in the leading ln(Q'^ / A^^^^) approximation outlined in Sec. 14.21 and found to be substantially 
enhanced as compared to the forward amplitude, 7*p 7*p [47 1. The DVCS cross section re- 
ported by the HERA experiments is in reasonable agreement with these predictions, as well as 
with the color dipole model of Ref. ll^ : see Ref. ll^ and references therein. The HERA data 
at small x are also well described by an NLO QCD analysis ll^ 1^ . in which the modeling 
of the input GPD's is a much more challenging problem than in LO, see Ref. fl65i for details. 
DVCS at small x and the closely related process of production of Z-bosons, j + p ^ Z + p, 
were also studied within the leading ln(xo/a;) approximation iH^ . 

To summarize, the HERA data on exclusive electroproduction of vector mesons clearly 
show the transition to the perturbative QCD regime for > 10 — 20 GeV^. This conclusion is 
consistent with the observation of color transparency phenomena in several other processes. It 
establishes the study of exclusive processes (x, and t-dependence of the cross section) as a 
way to extract detailed information about the interaction of small dipoles with hadrons, as well 
as about the generalized parton distribution in the nucleon. 



4.4 Transverse spatial distribution of gluons in the nucleon 

An important aspect of hard exclusive processes at small x is that they provide information 
about the transverse spatial distribution of gluons in the nucleon. It is contained in the Fourier 
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Figure 10: The average squared transverse radius of the gluon distribution in the nucleon, (p^) = 
/ <Pp p^Fg{x, p; Qeff)' ^ function of x, as extracted from J/ if) photoproduction data (Q^g = 
3 GeV^) at various energies. 

transform of the two-gluon form factor, Eq. dTSl) . 

F,{x,p-Ql^) ^ l^e^^^-^^F,{x,t = -Al;Ql^), (21) 

where p is a transverse coordinate variable. The function Fg(x, p; Qeg) is positive definite ll66l 
and describes the spatial distribution of gluons in the transverse plane, / d'^p Fg{x, p; Ql^) = 1. 

The convergence of the t-slopes of p and J/'ip production at large (see Fig.|9l) demon- 
strates that the t-dependence of the differential cross section is dominated by the two-gluon 
form factor. The two-gluon form factor can thus be extracted from the J/ 4' photoproduction 
data (Qlff ^ 3 GeV^), with small corrections (~ 10%) due to the finite transverse size of the 
J/ip meson. This process has been measured over a wide range of energies; see Refs. lIFTll^ 
for an overview of the data. At fixed-target energies, x ~ 10^^, the t-dependence of the data is 
well described by a two-gluon form factor of dipole form, 

Fg = (l-t/mD'^ ml = 1.1 GeV' (x ~ 10"^), (22) 

where the parameter, m^, is close to that in the dipole fit to the axial form factor of the nu- 
cleon. This corresponds to a narrow spatial distribution of gluons in the transverse plane, 
with an average transverse radius (p^) = ^ 0.28 fm^, see Fig. [lOl At HERA ener- 

gies, X ~ 10^^ — lO^'^, the average radius is larger, (p^) ^ 0.35 fm^. It also exhibits a slow 
growth with ln(l/a;), with a slope, a', significantly smaller than the value for soft interactions. 
The J/ip photoproduction data from HI give a;^^ = 0.08 ± 0.17 GeV"^ [69|, the ZEUS elec- 
troproduction data a{^^^ = 0.07 ± 0.05(stat)lo:o4(syst) GeV~^ JTOll, which should be compared 
to a' ~ 0.25 GeV^^ for pp elastic scattering. This reflects the suppression of Gribov diffusion 
for partons with large virtualities, see the discussion in Sec. 12.21 

The change of the nucleon's average transverse radius between x ~ 10^^ and 10^^ can 
naturally be explained by chiral dynamics. Pions in the nucleon wave function carry momentum 
fractions of the order m^/miv. For x > m.„/mM the pion cloud does not contribute to the gluon 
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distribution, and the two-gluon form factor is similar to the nucleon axial form factor, which 
also does not receive contributions from the pion cloud. For x <^ m^/mj^, the pion cloud 
contributes to the gluon distribution and leads to an increase of (p^) by 20 — 30%, see Fig.fTDI 

ED. 

It is worth emphasizing that for smaller x the increase of the transverse size should continue 
due to the Gribov diffusion. Indeed, a hard probe can interact with a parton of the soft ladder, 
responsible for the growth of the soft radius, if the soft parton's momentum fraction is suffi- 
ciently small compared to x. At very small x and fixed the rate of the growth should thus be 
comparable to that in the soft case ll72ll . No such effect is present in the BFKL model where the 
interaction of two small dipoles is considered. 

The change of the transverse spatial distribution of gluons in the nucleon with the scale, Q^f^, 
due to DGLAP evolution should generally be small [731 . For Qgff sufficiently large compared 
to the transverse spatial resolution, the parton decays happen essentially locally in transverse 
position. For fixed x, one finds that the transverse spatial distribution shrinks with increasing 
scale, because the distribution becomes sensitive to the input distribution (at the initial scale) at 
higher values of x, where it is concentrated at smaller transverse distances. 

4.5 Color transparency in hard processes with nuclei 

QCD predicts that the spatially small quark-gluon wave packets formed in hard 7*-induced 
scattering processes interact weakly with hadronic matter, because of the color neutrality of the 
photon. At sufficiently small x, where the cross section is proportional to the gluon density, 
cf. Eq. (El), one expects the ratio of the cross sections for 7* scattering from a nucleus and a 
single nucleon to be equal to the ratio of the respective gluon densities, a property known as 
generalized color transparency ifT^ OTl . Because with increasing gluon shadowing at fixed 
x disappears {cf. the discussion in Sec. l5.4l ). one further expects that 

a7ot^/(A(T7j^) ^1 {Q^ ^00; X fixed, small), (23) 

which is referred to as color transparency proper. Conversely, at fixed and decreasing x, 
the ratio in Eq. ( E^ should decrease owing to the more important role of nuclear shadowing, 
and color transparency phenomena should completely disappear at very small x, where QCD 
factorization breaks down. This is in contrast to the two-gluon exchange model of Refs. CH 
I75II . which neglects the space-time evolution of the dipole. In this model nuclear shadowing 
is obtained from exchanges of additional gluon between the current and target fragmentation 
regions, which is a higher-twist effect (oc l/Q"^) and disappears at large Q'^. 

The color transparency phenomenon has been directly observed in three experiments: 

• The total cross section for 7*^4 scattering increases with the atomic number as A" with 
a ~ 1, faster than the cross section for a hadronic projectile, see Ref. [i76il for a review of 
the experimental data. 
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• The cross section for coherent photoproduction of J/ ip mesons from nuclei increases with 
A much faster than that for coherent p meson production. The Fermilab E691 experiment 
11771 observed cr'y*+A^J/i'+A ^i.46 ^ 150 GeV. Color transparency predicts that 
the coherent cross section integrated over t is oc (A^/R^) ^ A'^^^. This A-dependence 
corresponds to the coherent sum of collisions from independent nucleons without absorp- 
tion. A somewhat faster A-dependence emerges because of the contribution of incoherent 
diffractive processes 

• The y4-dependence of the cross section for coherent dijet production in pion-nucleus 
scattering is anomalously large, as predicted in Ref. lfT6ll . The Fermilab E791 experiment 
IITS I observed a dependence oc A^-^^ at E^^ ^ 600 GeV, similar to that in coherent J/ip 
production. Note that the conventional Glauber approximation predicts A^^^. Further- 
more, the observed dependence of the cross section on the pion momentum fraction and 
the jet transverse momentum is well consistent with the perturbative QCD prediction of 
Ref. llT9l . Notwithstanding the fact that the absolute cross section has not been measured, 
this is probably the first experimental observation of the high-momentum tail of the pion 
wave function as due to one-gluon exchange. 

To summarize, there exists strong experimental evidence for color transparency in high-energy 
scattering. This phenomenon could be the basis for new "non-destructive" methods of investi- 
gating the microscopic structure of hadrons and nuclei in the future. 

5 Diffraction in Yp scattering 

5.1 QCD factorization for hard diffractive processes 

Measurements of DIS at HERA have established the existence of a class of events in which the 
proton is observed in the final state, with a small invariant momentum transfer, t, and a hadronic 
system of invariant mass Mx -C W'^ is produced with a rapidity gap relative to the proton. In 
a frame in which the nucleon is fast-moving (i.e., in parton model kinematics) such processes 
are characterized by the fractional energy loss of the proton, xp = (E^ — E,p/Ep, and the 
transverse momentum transfer, A^, with t = — (A^^ + x'jpm'j^) /(I — xp). In analogy with the 
corresponding phenomenon in hadronic collisions one refers to such processes as diffractive, 
although a priori the dynamics is not governed by soft physics. 

Following suggestions of earlier works, a formal QCD factorization theorem was proved in 
Ref. |I79l|80J for hard processes of the type 

7* + p ^ h + (rapidity gap) + X, (24) 

where X is either an inclusive state, or a state with extra hard activity (dijet production, heavy 
quark production, etc.), see Fig. [TTl Similar to inclusive DIS, processes (0^ with a given 
hadron h in the target fragmentation region are characterized by so-called conditional parton 
distribution functions, Q^, Xh, t), which are independent of the hard process and satisfy 
the same DGLAP evolution equations for fixed Xh and t. Here {3 = x/ {1—Xh) = I (Q^+^x) 
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Figure 11: Factorization in diffractive DIS. The amplitude in the dashed blob, multiplied with 
its complex conjugate and summed over final states, defines the diffractive PDF The internal 
structure of the dashed blob illustrates the assumption (l25t . which does not follow from the 
QCD factorization theorem. 

is the fraction of the light-cone momentum of the target available for hard interactions in hadron 
h. Most of the current studies focus on diffractive kinematics, where h = p and 1 — Xh = xp < 
0.01. The conditional parton distribution functions in this case are referred to as diffractive 
parton distribution functions (dPDF's), and denoted as /j^. In current data analysis it is usually 
assumed that the dependence of the dPDF's on xp, t and /5, can be factorized as lISTIl 

fl'iP^Q^xjp^t) = fjp/pixjp,t) ff{P,Q'), fipUxjp,t) = /(t) x^'"-W+2_ ^25) 

This assumption is inspired by the soft Pomeron exchange model (which does not follow from 
the QCD factorization theorem) and referred to as Regge factorization, see Fig. [TTl An addi- 
tional term can be added to Eq. (E51) in analogy to non- vacuum exchange in soft physics; it gives 
a small contribution below xjp ~ 0.01 and dominates at xp > 0.05. 

Extensive studies of hard diffractive channels have been performed at HERA. The inclusive 
diffractive cross section was measured both integrated over t (the so-called diffractive structure 
function, F^^^^), and as a function of t for a limited range of xjp. These measurements are 
mostly sensitive to the quark dPDF, while the gluon dPDF enters through scaling violations. 
Diffractive dijet production for real and virtual photons, as well as diffractive charm production, 
primarily probe the gluon dPDF The analysis of these data on the basis of QCD evolution 
equations has led to the following conclusions: 

• The data at > 4GeV^ are described by the universal dPDFs, consistent with the 
factorization theorem. 

• fgiP,Qo) ^ J2q fqiP^Qo) for the studied range of /9. This conclusion was initially 
based on the weak scaling violation for F^^^\(3, Q"^) for large (3, and was later confirmed 
by the studies of diffractive dijet production and charm production. However, the latter 
processes have so far been treated only in the LO approximation, and one should await 
the NLO analysis before drawing final conclusions. 
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Figure 12: Probability of diffractive scattering from anti-n quarks (left) and gluons (right), 
extracted from a fit to the HI data, see Ref. lIH^ (updated by V. Guzey). 

• The data are consistent with Regge factorization, Eq. (E51) . although the xp dependence 
is faster than in soft physics. The analysis of ZEUS and HI diffractive data finds ap^t = 
0) = 1.2 ± 0.07 and increasing with Q"^, which should be compared with the energy 
dependence expected in soft hadronic collisions, ajp ^ 1.1. 

• The absolute probability of diffraction in 7* +p scattering is of the order of 10% for mod- 
erate Q^, and thus of the same magnitude as in soft pion-nucleon collisions. However, 
the rate of increase of the diffractive cross section with energy for fixed Mx and is 
significantly faster than that of the total DIS cross section. 

Another interesting characteristic of diffractive DIS is the probability of diffractive scatter- 
ing depending on the type of parton coupling to the hard probe fT2\ . 

Pjix, Q2) = jdtj dxjp ff{x/xjp, Q\ xjp, t) I f,{x, Q2). (26) 

This ratio cannot exceed the value 0.5, which corresponds to the unitarity (black disk) limit 
(BDL), cf. the discussion in Sec.0 Using the HI fit to the diffractive DIS data (see Fig.fT^. 
we find Pg > Pg, and Pg{x ~ lO^^) ^ 0.4 (0.3) for = A (10) GeV^. That is, quark induced 
diffraction is small, whereas gluon induced diffraction is close to the maximum value allowed 
by unitarity. We shall return to this point in our discussion of the profile function for the dipole- 
nucleon interaction in Sec. |6l Note that the HI fit is based on the data at x > 10~^. The fact 
that it leads to Pg > 0.5 at smaller x indicates that the HI parameterization should break down 
near the upper end of the HERA energy range. 

5.2 Space-time picture of hard diffractive processes 

To understand the observed pattern of hard diffraction, it is instructive to consider the space- 
time evolution of such processes in the target rest frame. Such studies reveal new information 
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about the interaction of small-size qq as well as qqg . . . g configurations with hadronic matter. 
In particular, the ratio of the diffractive to the total cross section probes the interactions of such 
configurations with the target without reference to the probe which created them. 

An immediate consequence of the QCD factorization theorem is that, in the target rest frame, 
the number of components in the photon wave function evolves with Q^. While at low an 
approximation by a few components in the photon wave function may be reasonable, it is defi- 
nitely inappropriate for large Q'^. This point is illustrated by the following example: Diffractive 
processes induced by longitudinally polarized photons are a leading-twist effect. If, however, 
all but the qq component of the photon wave function were neglected, one would erroneously 
conclude that diffraction is a higher-twist effect in this channel, because the transverse size of 
the longitudinal photon is d'^ oc l/Q^. The proper dependence is restored by the qqg . . . g 
configurations in the photon wave function. 

At low Q^, aligned jet model-type configurations of large transverse size dominate in the 
wave function of the projectile photon, cf. Sec. 12.21 Such configurations interact with a hadronic 
cross section, ~ (Jtot(7rA^), and thus have a significant probability to rescatter elastically. If Q'^ 
increases, these configurations cannot be effectively produced without emission of gluons. Be- 
cause these gluons are predominantly emitted coUinearly, they do not change the transverse size 
of the diffracting system, and hence the probability of elastic rescattering. Assuming smooth 
matching between the strength of interaction in the perturbative and nonperturbative regimes, 
models reasonably describe the data on hard diffraction in ep scattering if Pg > Pg at the initial 
scale, Ql ll83l [84ll . Configurations of size cP > \/Ql, where A ~ 9 (c/ the discussion in 
Sec. 13.21) should be included in the definition of the dPDF at the initial scale. For Ql = A GeV^, 
this includes rather small transverse sizes, for which the cross section increases with energy sig- 
nificantly faster than at the soft scale, which is consistent with the trend of the HERA data men- 
tioned in Sec. 15. 11 Also, because the qqg configurations have masses considerably larger than qq 
configurations, they should manifest themselves in diffraction at relatively small (3 < 0.5. Indi- 
cations of diffraction into non-aligned jet final states were indeed found in a number of HERA 
experiments [85|. To summarize, it appears that hard diffraction at HERA with ~ 4GeV^ 
represents the border between the high-Q^ region where leading-twist QCD dominates and the 
low-Q^ region where higher-twist effects become important. 

A quantitative analysis of the HERA data within the gluon dipole picture indicates that the 
interaction of gluon dipoles at HERA energies is rather close to the BDL cf. the discussion 
in Sec. |5| However, due to our inability to build an effective trigger on the interaction of gluon 
dipoles with d > 0.3fm it is difficult to observe this effect directly in the experiments. 

To conclude this discussion, we briefly want to comment on the assumption of Regge fac- 
torization, Eq. (E51 . Because strong deviations of the energy dependence of diffraction from the 
soft regime are observed, there is a priori no reason for the validity of this assumption. Several 
effects are likely to contribute to the breakdown of Regge factorization: (a) different energy 
dependence of the cross section for diffraction of configurations of different transverse size, 
(b) emission of gluons by qq dipoles, which at smaller x occurs at large coherent lengths {cf. 
the discussion in Sec. 12. 2I ). and (c) soft screening effects, which become more important with 
increasing energy (these effects were observed in soft diffraction ll86ll87l '). and which should 
be different for the various diffractive configurations, as they interact with different strengths. 
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One way to probe the degree of validity of Regge factorization is to check that the hadrons 
produced in the photon fragmentation region do not "talk" with hadrons in the target fragmen- 
tation region, i.e., that there are no long-range correlations in rapidity. Such a factorization 
ignores the existence of color fluctuations, which lead to processes in which the proton also 
dissociates. If Regge factorization were valid, the probability of dissociation would not depend 
on the properties of the state into which the virtual photon has diffracted. However, if screening 
effects are present and non-universal (for example, due to the different strength of interaction in 
quark and gluon induced diffraction), Regge factorization should be broken. 



5.3 Diffraction and leading-twist nuclear shadowing 

A direct relation between diffraction in high-energy hadron-hadron collisions and the nuclear 
shadowing effect in hadron-nucleus collisions was derived by V. Gribov lISSllS^ . in the approx- 
imation where the nucleon radius is considerably smaller than the mean intemucleon distance 
in nuclei. The same reasoning in conjunction with the leading-twist approximation for hard 
diffractive processes allows one to calculate nuclear shadowing of PDF's in light nuclei [i72J . 

Application of the Abramovsky-Gribov-Kancheli cutting rules f90l, which relate the shad- 
owing phenomenon for the total cross section of high energy processes to that of partial cross 
sections (such as for diffraction, multiparticle production, etc.) and are valid in perturbative 
QCD (see Ref. 11911 for a recent discussion), explicitly demonstrates that the interference of 
the amplitudes of diffraction from a proton and a neutron leads to a decrease of the total cross 
section for 7*D scattering. When combined with the factorization theorem for inclusive diffrac- 
tion, one can calculate the modification of nuclear PDF's at low values of Bjorken x [72|. In the 
limit of low nuclear thickness, the nuclear shadowing corrections to the nuclear parton densities 
are given by 

-T = Jj/N{X,Q > " 2 COS[X]pmN[Zi - Z2)\ 

X ^^^2 fflN{l3^Q'^^^lP^'>)\k^=Q pAib,zi) pAib,Z2), (27) 

where fj/Nix, Q"^) is the usual parton density in the proton, f^jj^if^, Q^, xjp, t) the diffractive 
parton density (see Sec. 15. II) . and pAi^) is the nucleon density in the nucleus with atomic num- 
ber A. The momentum transfer, t, is given by — t = {k^ + {xjpmpf)'^) / (1 — X]p), where k± 
is the transverse component of the momentum, transferred to the struck nucleon, and /3 = 



x/xip. In Eq. the factor (1 - r/2)/(l + r]'^), where r] = -tt /2d\Yi{^ f^^) / d\n{l/x]p) = 

71 /2 [a]p{t = 0) — 1], accounts for the real part of the amplitude of diffractive scattering [92 1. 
One can easily modify Eq. (l27t to include the dependence of the diffractive amplitude on t. 
Obviously, both the left- and right-hand side of Eq. (ITTt satisfy the QCD evolution equations 
in all orders in a^, and this relation does not depend on the renormalization scheme. These 
expressions represent the model-independent result for leading-twist nuclear shadowing in the 
low-density limit. 

In the case of heavy nuclei one may with good accuracy neglect the fluctuations of the 
strength of interaction in the hadron component of the photon wave function. This approxi- 
mation makes it possible to extend the above formulas to the case of heavy nuclei [|72j, .93J . 
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Numerical studies of shadowing using Eq. dTTl) and the corresponding expression for the total 
DIS cross section for heavy nuclei found large leading-twist shadowing effects for quark and 
gluon distributions, with gluon shadowing being larger up to rather high values of . The 
latter effect can be traced to the higher probability of gluon-induced diffraction as compared to 
quarks, see Fig.fT2l 

The connection between diffraction and nuclear shadowing does not depend on the twist 
decomposition of the cross section, and was successfully applied also to data on nuclear shad- 
owing of F2A at intermediate Q^, see e.g. Refs. Il94ll95l . One can use this to estimate the relative 
importance of leading-twist and higher-twist nuclear shadowing at < 2 GeV^, using exper- 
imental information on the leading-twist contribution to the diffractive cross section at these 
values of . One finds that a significant higher-twist contribution to diffractive DIS originates 
from p meson production. This implies that a significant fraction (~ 40%) of the nuclear shad- 
owing observed in the experiments at CERN and Fermilab (see Ref. 1761 for a review) may be 
due to higher-twist effects W^ . 

Recently, leading-twist nuclear shadowing was included in the initial conditions for the 
small-x evolution in the McLerran-Vegnugopalan model ll96l . A distinctive feature of this 
model is that the invariant masses in the nuclear vertex of the BFKL ladder should be very large 
as compared to . That is, small /5 should dominate in the integral, in analogy to Eq. (l27t . A 
numerical analysis of gluon shadowing using the HERA dPDF's finds that the region /3 < 0.1 
becomes important only for x < 10"''. Thus, the assumption of the dominance of large diffrac- 
tive masses may give rise to important dynamical effects at the next generation of accelerators. 

5.4 Implications of nuclear shadowing for heavy-ion collisions 

The typical x-values relevant for semihard production of hadrons in heavy-ion collisions de- 
creases with energy as xa ~ 2p_l/ ^/snn for central rapidities, and much faster, oc 1/s, for 
the fragmentation regions (sat at is the squared center-of-mass energy of the effective nucleon- 
nucleon collisions). For central rapidities and p± > 2GeV, gluon shadowing is still a small 
correction at RHIC. However, it will be a large effect at LHC for a wide range of p±, because 
the relevant xa are much smaller than 0.01 . The expected suppression of jet production is given 
by the factor [Ga{xa, P±)/^G{xa, P±)]^, where Ga and G are the gluon densities in the nucleus 
and the nucleon, respectively. This factor can be of the order of 1/4 II82I . 

Because the current RHIC detectors have rather limited forward coverage, they have lim- 
ited sensitivity to small-x phenomena. One exception is J/ip production, which, if interpreted 
within perturbative QCD, probes x down to 0.003. The observed suppression of the J/'tp yield 
is consistent with the estimates of Ref.[93|, see Ref. [97 1 for a review. The A-dependence of 
forward high-p^ hadron production was studied by BRAHMS [98 1. Although at large rapidities 
small X contribute to the high p± spectra, the QCD analysis indicates that average s are ~ 0.03 
[1991 . Consequently, the yields are practically not sensitive to the shadowing effects, or, more 
generally, to any initial-state modifications of the nucleus wave function. Final-state interac- 
tion effects which could explain the data are nonperturbative contributions to the production of 
leading hadrons, due to coalescence of spectator partons and the relatively small energy losses 
in the initial and final state (on the scale of 3% ). 
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6 Black-disk limit in dipole-hadron interactions 
6.1 Violation of the leading-twist approximation at small x 

A fundamentally new dynamical effect expected at high energies is the unitarity limit, or black- 
disk limit (BDL), in the interaction of a small dipole with hadronic matter. We now describe 
and quantify this effect, using the information gathered in our studies of inclusive, exclusive 
and diffractive DIS in Sections QUI and 13 

A simple argument shows that the twist expansion for the cross sections of hard processes 
breaks down at sufficiently small x. QCD factorization predicts that the total cross section for 
DIS at fixed Q"^ increases with decreasing x as (Xtot oc xG{x)/Q'^. At the same time, the cross 
section of elastic dipole-hadron scattering (which corresponds to the production of diffractive 
states with masses Mx oc 1/d oc Q) grows much faster, adiff oc (xCy /(BQ^), where B is the 
t-slope of the corresponding differential cross section [21J. Clearly, at sufficiently small x there 
is a contradiction — the total cross section should always be larger than that for any particular 
channel ll44ll . The resolution of this paradox is that the decomposition of hard amplitudes in 
powers of l/Q^ becomes inapplicable at sufficiently small x. 

In order to quantify the onset of the new regime, it is instructive to consider the effects 
of unitarity (conservation of probability) on a purely theoretical scattering process, namely 
the scattering of a qq (quark-antiquark) or gg (gluon-gluon) dipole of small transverse size, 
d, from a hadronic target. Neglecting other constituents in the dipole is justified in a wide 
kinematic range by the smallness of the coupling constant; large terms oc asln(xo/x) arise 
only from interactions at large rapidity intervals. The invariant amplitude for dipole-proton 
elastic scattering is a function of the invariants s = W"^, and t. We write is as a Fourier integral 
over the dipole-proton impact parameter, b. 



where r*(s, b) is the so-called profile function. Making use of unitarity, one can express the 
total, elastic, and inelastic (total minus elastic) cross sections in terms of the profile function as 



In the situation where elastic scattering is the "shadow" of inelastic scattering, the profile func- 
tion at a given b is restricted to values iV^^^s, b)\ < 1. The value r*(s, b) = 1 corresponds to 
complete absorption at a given impact parameter, the so-called black disk limit (BDL). 

The proximity of r*(s, b) to unity is an important measure of the strength of the interaction 
of the dipole with the proton. As outlined in Sec. 12.21 the analysis of 'j*N scattering in the target 

'^In non-relativistic quantum mechanics, the scattering of a particle from a completely absorptive sphere is 
referred to as the "black-body limit." In contrast, the high-energy limit of scattering amplitudes in QCD is essen- 
tially two-dimensional, with the radius of interaction increasing with energy. It is thus natural to refer to the limit 
of complete absorption as the "black disk limit." 





(28) 




(29) 
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d = .1 (fm) d = .3 (fm) d = .5 (fm) 



Figure 13: The profile function of dipole-nucleon scattering, r"'^, as a function of the impact pa- 
rameter, b, for various values of the dipole size, d, and x, as obtained from a phenomenological 
estimate (see text). 

rest frame allows one to determine with reasonable accuracy (LO approximation) the total cross 
section for the scattering of a qq dipole from the proton. Combining this with information 
on the transverse spatial distribution of gluons in the nucleon, obtained from measurements 
of the t-dependence of exclusive J/ photoproduction and other hard exclusive processes (cf. 
Sec. 14.41) . we can calculate the profile function for dipole-nucleon scattering |73|. A sample of 
the results of Ref. lllQOII is presented in Fig. [131 Here, x = Q'^/s, where Q'^ = X/<P ^ 9/(P 
is the characteristic virtuality corresponding to the dipole size. One sees that up to the top 
of the HERA energy range, x ~ 10^'', the profile function for a qq dipole remains small for 
dipole sizes d < 0.3 fm, corresponding to > 4GeV^, the value usually used as a starting 
point for DGLAP evolution 11011 1. Even for larger dipole sizes, d ~ 0.5 fm, the fraction of the 
cross section due to scattering with ReF > 0.5 remains small. This shows that the BDL is not 
reached in inclusive DIS at HERA energies, in agreement with what we observed in Sec. 13.11 
and l3.2[ However, the DGLAP evolution starting from a qq dipole generates gg dipoles, whose 
interaction at leading twist is larger by a factor of 9/4, cf. Eq. ©, and thus approaches the BDL 
much earlier. The theoretical estimate shown in Fig.[T3]indicated that the strength of interaction 
of gg dipoles is close to maximal at the top of the HERA energy range, for a wide range of b, and 
dipole sizes corresponding to < A GeV^. This implies that in gluon-induced interactions at 
HERA at such Q"^ the probability of diffraction should be close to 1/2 — exactly as we found 
in the analysis of HERA diffractive data in Sec. 15. 1[ 

To summarize, one may expect that the leading-twist approximation for DIS breaks down 
for < 4GeV^, especially in the gluon sector. Unfortunately, there are no readily available 
probes of gluons at low Q'^, except possibly the longitudinal cross section, a^. Without such 
measurements, it is impossible to determine whether the successful DGLAP fits to the HERA 
data down to ~ 1 GeV^ are an artifact of using an essentially arbitrary gluon density at low 
Q"^ and low x; this function is practically not constrained by the data at larger Q'^, where it is 
dominated by DGLAP evolution from larger x. Also, it is worth emphasizing that although 
F'^*'(6) for gluon dipoles (and, at higher energies, also for quark dipoles) reaches values close to 



34 



unity, the actual deviation of the nucleon structure function from the DGLAP fits may still be 
rather small, as the contributions from scattering at large b remain dominant. 

6.2 Theoretical issues in describing the black-disk limit 

The analysis of the full DGLAP evolution equation and the resummation approaches shows that 
the cross section for the scattering of a small dipole increases effectively as x^", with n > 0.2 
for Q"^ ~ few GeV^, and somewhat faster at larger Q^. [At very large Q'^ and extremely high 
energies — well above the LHC range — the resummation approaches predict that n oc asiQ"^), 
and thus n 0.] Within these approximations, the t-slope of the differential cross section for 
elastic dipole-nucleon scattering, B, increases with energy rather slowly, cf. the discussion in 
Sec. 14.41 Thus, the cross section for dipole-nucleon elastic scattering increases with energy 
faster than the total cross section, and unitarity is violated for this hard processes within the 
leading-twist approximation (see the discussion in Sec. 16. II) . Probability conservation will be 
violated at fixed impact parameter, b, in a region corresponding to a disk of finite size. At 
large enough b, the interaction is too small to violate the leading-twist approximation. Note 
that in soft interactions the increase of the elastic cross section does not necessarily lead to a 
contradiction with unitarity, because soft interactions generate also a', and therefore give rise 
to an increase of the t-slope as 5 oc In s. 

The conventional assumption is that, beyond the leading-twist approximation, taming of 
the growth of cross sections occurs due to the shadowing phenomenon (this follows e.g. from 
V. Gribov's reggeon calculus). Specific to this phenomenon is that bare particles may expe- 
rience only one inelastic collision, but any number of elastic interactions, without changing 
trajectory. The behavior of the amplitudes for high-energy processes in QCD differs from that 
given by the eikonal approximation in non-relativistic quantum mechanics because of the ne- 
cessity to account for the non-conservation of the number of bare particles. Application of the 
Abramovsky-Gribov-Kancheli cutting rules |90| shows that the requirement of positive proba- 
bilities for total cross sections, single particle densities, etc., impose serious restrictions on the 
dynamics in the case of cross sections increasing with energy. To satisfy these requirements in 
a series of multiple rescatterings in which consecutive terms have alternating signs, the effec- 
tive number of constituents in the dipole wave function should increase with energy. Thus, the 
increase of the cross section with energy leads to resolution of constituents in the wave func- 
tions of the colliding particles, and therefore to an evolution of final states. Evolution and gluon 
emission by dipoles are the key for generating multiple inelastic collisions without violation 
of causality and energy-momentum conservation. The evolution of a dipole in time manifests 
itself in the expansion of the system, emission of gluons, transitions between components con- 
taining different numbers of bare particles, change of impact parameters in the intermediate 
states, and the related effect of the cross section for inelastic diffraction exceeding the elastic 
cross section for the scattering of small dipoles, cf. the discussion in Sec. |5l In this regime, 
the concept of a parton density of the target cannot be defined in a model-independent way, 
because the parton distributions in the dipole and the target are intertwined and not restricted by 
probability conservation. 

At energies where the dipole cross section becomes comparable or even larger than 27r_R^, 
the whole picture of rescattering becomes inconsistent if the radius squared of the interaction 
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is not proportional to the dipole cross section. For example, in this case the Glauber approxi- 
mation for hadron-deuteron scattering would lead to negative total cross sections. Fortunately, 
in QCD the wave function of a fast projectile contains many partons. This fact, combined with 
the increase of the dipole cross sections with energy, is sufficient to ensure complete absorp- 
tion for central collisions, without detailed knowledge of the hadron wave function [102|. To 
illustrate the rapid onset of complete absorption for central collisions related to the increase of 
the number of constituents, we adopt here a simple approximation, namely that the perturbative 
QCD (pQCD) description of dipole-proton scattering works up to Re r*(6, x) = 1, and that 
r*(6, x) = 1 if the pQCD formulas lead to Re r*(6, x) > 1. That is. 



Many of the models currently discussed in the literature use the elastic eikonal approxima- 
tion to describe the taming of the increase of the dipole-hadron cross section with energy as due 
to the shadowing phenomenon, see e.g. Refs. II103LI1041I for a nucleon target and Refs. 1 105 1 for 
a heavy nuclear target. These models assume that taming becomes significant for r*(6) > 0.5, 
i.e. at significantly larger x than the values where unitarity is explicitly violated in the pQCD 
approximation. Early taming results in a slow approach to the unitarity limit, T^P{b) = 1. Obvi- 
ously, these conclusions are model dependent, as such models neglect most of the QCD effects 
mentioned above. 

The condition of the BDL for dipole-hadron scattering, r*(6) = 1, expresses the complete 
"loss of memory" of the cross section on the structure of the projectile and the target, and of 
the value of the running coupling constant, in a finite region of transverse space. It reflects 
the breakdown of two-dimensional conformal invariance (which is the basis of approximate 
Bjorken scaling in DIS) because of the appearance of a dynamical scale related to the high 
gluon density and the radius of the transverse distribution of gluons. The qualitative departure 
from pQCD dynamics in the BDL cannot be explained as a soft interaction effect. This can be 
understood when considering collisions of two small dipoles of same size near the BDL, e.g. 
7*((5^)-7*(Q^) scattering at sufficiently large Q^, in which soft interaction effects are under 
control and negligible. An interesting question is whether, from a general perspective, the 
appearance of this new scale corresponds to a spontaneous breakdown of conformal symmetry, 
or is related to the conformal anomaly. 

6.3 High-energy limit of nuclear and hadronic structure functions 

The approach to the BDL in dipole-hadron scattering at high energies has interesting implica- 
tions for the theoretical behavior of hadron and nuclear structure functions at extremely high 
energies, which is subject to the Froissart bound. 

We consider the scattering of a virtual photon from a heavy nucleus (radius Ra) at high 
energies as a superposition of the scattering of dipoles of different sizes. The interaction at 
impact parameters b < Ra will be black for dipoles with sizes larger than some critical size, 
d > d(x), leading to a contribution to the cross section [cf. Eq. dTTTl l 



Rer*(6,x) 



= Rer*(6,a;)pQCDe[l-Rer'^^(6 
+ e[Rer'^^'(6,x)pQCD-l]. 



(30) 



cr 




(31) 
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Because the profile function of dipole-nucleus scattering increases like a power of energy in the 
region where it is < 1, one concludes that d{x) oc x™, with m > 0. When the nuclear radius 
significantly exceeds the essential impact parameters in 7*A^ collisions, one has 

F2aM') = ^(1:4) (2>rfli)ln^. (32) 

where xq{Q'^) slowly decreases with increasing Q'^.^ For illustrative purposes, we neglect here 
the contributions from peripheral collisions, which grow with the atomic number as A^^^. Al- 
though formally these contributions increase with energy faster than those from central colli- 
sions (which are nearly energy-independent), they are still comparatively small at all achiev- 
able energies. The gross violation of Bjorken scaling in Eq. (B^ . F2A oc Q"^ ln(xo/x), and 
the numerical coefficient follow from the normalization of the photon wave function to the Q^- 
derivative of the photon polarization operator, as opposed to unity as for hadron wave functions. 
Equation (l32t represents a QCD modification of the Gribov BDL [89 1, which assumed that all 
configurations in the virtual photon with masses < 2mNx/RA interact with the heavy 
nucleus with maximum strength. 

In DIS from a proton target scattering at large impact parameters is always important in 
the regime where the pQCD interaction becomes strong. Indeed, based on the studies of the 
transverse spatial distribution of gluons in hard exclusive processes (see Sec. 14.41) one expects 
that r*(s, b) oc exp(— /i6) at large b, with /i ^ nig for moderately small x [cf. Eq. (E^ l. and 
/i 2m^ in the limit of infinitely large energies. It follows from the unitarity bound, Eq. (E^ . 
that that the essential impact parameters increase with energy as 6^ oc In^ [(t(s, d) / (SnB)], 
where B is the t-slope of the differential cross section of dipole-nucleon scattering, which 
is almost energy-independent within the leading-twist approximation. The cross section of 
dipole-nucleon scattering therefore increases with energy as [1IO6I 



-2 1 2 cr{s,d) 

a (X u m . (33) 

In general, this behavior differs from the Froissart limit for soft hadronic interactions, because 
of the more complicated dependence of the dipole cross section, a{s,d), on the energy, as 
described by the resummation approaches. To simplify the formulas, below we shall use the 
observation that effectively a(s, d) oc s"("'), with n > 0.2 for small d. This approximation leads 
to the limiting behavior familiar from soft hadronic interactions 1 3 1 . The leading asymptotic 
term in x for fixed for the nucleon structure function is 

F^M') = 7^l>:4kln- « (34) 

where two logarithms originate from the dipole-nucleon cross section, and one from the integral 
over the photon wave function, similar to the case of scattering from nuclei. 




^Since in the BLD the masses of the intermediate states are much larger than Q^, the coherence length is much 
smaller than the naive estimate, ^coh ~ l/{2mNx). If the gluon density in the approach to the BDL grows as 



X 



^HQ )^ one expects that Icoh oc l/{mpfX^ ^ 
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In hard exclusive processes in scattering, the approach to the BDL implies that the t- 
slope increases with energy oc ln^(s/so) {cf. also Sec. l7.4l ). A promising strategy in searching 
for BDL effects would be to extract partial waves for small impact parameters from the cross 
sections of processes such as DVCS from the nucleon, p-meson production II107L coherent pho- 
toproduction of high dijets from the nucleon, and coherent photoproduction of J /if), -ip', tp" 
mesons. This would allow one to probe small b, where the gluon density is maximal and unitar- 
ity effects should manifest themselves early. Another possible strategy is to study the structure 
functions of heavy nuclei, where due to higher gluon densities the BDL effects are enhanced by 
a factor A^^^l 108 1. Note that this enhancement is partly compensated by nuclear shadowing; cf. 
the discussion in Sec. 15.31 

It is interesting also to explore the behavior of nuclear structure functions at extremely high 
energies, where the radius of the 7*A^ interaction becomes comparable to or even exceeds the 
nuclear radius. In this case, first the edge of the nucleus contributes terms oc A^^^ ln(xo/x)^ 
to the cross section. Ultimately, for s — > cx) and fixed Q^, one would reach the universality 
regime where F2yi(x, Q'^)/F2p{x, Q"^) 1 II102II . However, the relevant scale is comparable to 
the gravitational scale. 



6.4 Black-disk limit in hard diffractive scattering from heavy nuclei 

An important consequence of the BDL is that, in diffractive scattering at high energies, non- 
diagonal transitions between diffractive eigenstates are forbidden |89|. This follows from the 
orthogonality of the eigenstates — if every configuration in the projectile interacts with the 
same strength, the relative proportion between different configurations remains the same. This 
implies that half of the nuclear DIS cross section should be due to diffraction, with the nucleus 
remaining intact, and a "jetty" diffractive final state resembling that of e^e~ — > hadrons. In 
contrast, in the leading-twist approximation this cross section should be negligible. 

At the onset of the BDL regime, where the contributions from configurations in the virtual 
photon interacting with the BDL strength and those for which pQCD is applicable are compa- 
rable, one can calculate the differential cross section for diffraction to final states of small mass, 
Mx, for which the interaction is already black, 

dFf'^''{x,Q\Ml) ^ ttRI Q^Ml da{e+e- ^ X)/dQx 

dMld^x 127r3 {Ml + Q2)2 ^{e+e- n+fi-) ' ^ ^ 

This shows a much slower decrease with than in the leading-twist approximation, and cor- 
responds to "jetty" final states (mostly diffraction to qq and qqg jet states). The earliest signal 
for the change of dependence should be in p meson production, where the Q'^ dependence 
of the dominant longitudinal cross section should change from (see Sec. 14.31) to l/Q^. 

Theoretical studies show that at HERA kinematics the fraction of the cross section due to 
diffraction should be much closer to 1/2 for nuclear targets than for the proton 1 109 1. The use 
of nuclear beams would greatly facilitate the exploration of the BDL regime. Possibilities for 
such measurements in ultraperipheral collisions at LHC will be discussed in Sec.|5| 
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7 Small-x dynamics in hadron-hadron collisions 



7.1 Transverse radius of hard and soft interactions 

Our studies of 7*p scattering at HERA have taught us several important lessons about the gluon 
density at small x, the transverse spatial distribution of gluons, and about the interaction of 
small-size color singlets with hadrons (see Sec.[l]for a summary). We now explore the impli- 
cations for the physics of pp/pp and pA collisions at high energies. 

In hadron-hadron collisions, hard processes arise from binary collisions of partons in the 
colliding hadrons, in which a system of large invariant mass, ^ ^qcd' produced. Ex- 
amples are the production of dijets, dilepton pairs (Drell-Yan process), and the production of 
heavy particles such as Higgs bosons or SUSY particles. Such hard partonic processes are gen- 
erally accompanied by a rich spectrum of soft interactions, which dominate the total hadronic 
cross section and determine the overall characteristics of hadron production in the final state. 
Understanding the interplay of hard and soft interactions is the main challenge in describing 
hadron-hadron collisions with hard processes. 

A crucial observation in studies of hard exclusive processes in 7*p scattering is that gluons 
with significant momentum fraction (x > 10^^) are concentrated in a transverse area much 
smaller than that associated with the nucleon in pp elastic scattering at high energies, which 
is dominated by soft interactions. The difference between the two areas becomes more pro- 
nounced with increasing energy. When considering the production of a system of fixed mass, 
M^, in collision of partons with xiX2 = M'^/s, the transverse area of the hard partons grows 
with energy as (p^) = a[^^^ In s, whereas the transverse area for soft interactions grows at a 
much faster rate, a'^^f^ ^ Q^Lrd- The cause of this difference is the suppression of Gribov diffu- 
sion for partons of large virtuality, as described in Sec. l2.2[ Thus, in high-energy pp collisions 
one is dealing with an "onion-like" transverse structure of the nucleon (two-scale picture), see 
Fig. nil. 

The two-scale picture of the transverse structure of the nucleon implies a classification of 
pp/pp events in "central" and "generic" collisions, depending on whether the transverse areas 
occupied by the large-x partons in the two protons overlap or not, see Fig.fTIb |73|. Generic 
collisions give the dominant contribution to the overall inelastic cross section. Hard processes, 
such as heavy particle production at central rapidities, will practically happen in central colli- 
sions only. (Obviously, in these collisions multiparton interactions due to the small-x gluon 
fields are strongly enhanced, giving rise to the dynamical effects described in the following 
subsections.) 

To quantify the distinction between generic and central collisions, we estimate the distribu- 
tion of the probability for both types of events over the impact parameter of the pp collision, 
b. For generic collisions, the distribution is determined by the 6-dependent probability of in- 
elastic interaction, obtained via unitarity from the elastic pp amplitude in the impact parameter 
representation, Eq. (l29l) . We define a normalized 6-distribution as 

PUs,b) = [2RcTffis,b)-\TPPis,b)\'']/aUs), (36) 

where crin(s) is the inelastic cross section, which is given by the integral J cPb of the expression 
in the numerator. For collisions with a hard process, on the other hand, the 6-distribution 
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Figure 14: (a) The two-scale picture of the transverse structure of the nucleon in high-energy 
collisions, (b) The resulting classification of pp/pp events in "generic" and "central" collisions. 



is determined by the overlap integral of the distribution of hard partons in the two colliding 
protons (see Fig.fTib). 

P2{h) = j d^pi j d^P2 5^^\b- pi + P2) Fg{xi,pi)Fg{x2,p2). (37) 

Numerical estimates can be performed with our parametrization of the transverse spatial dis- 
tribution of hard gluons, see Sec. 14.41 which takes into account the change of the distribution 
with X and the scale of the hard process. The two ^-distributions are compared in Fig. fTSlfor 
Tevatron and LHC energies. For the hard process we have taken the production of a dijet with 
transverse momentum q±_ = 25 GeV at rapidity y = in the center-of-mass frame; in the case 
of Higgs production at LHC the distribution P2{b) would be even narrower. The results clearly 
show that events with hard processes have a much narrower impact parameter distribution than 
generic inelastic events. 

One expects that at LHC energies the rate of production of two pairs of jets will be very 
high. It is interesting to consider the ^-distribution also for the production of two dijets in two 
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Figure 15: Solid lines: Impact parameter distributions in events with hard dijet production, 
P2{b), in pp collisions at LHC (left panel) and pp collisions at the Tevatron (right panel), for 
a back-to-back dijet at zero rapidity with transverse momentum = 25 GeV. Long-dashed 
lines: Same for double dijet events, Pi{h). Short-dashed lines: Impact parameter distributions 
in generic inelastic collisions, obtained from the parameterization of the elastic pp amplitude 
of Islam et al. [JIOJ . The plots show the "radial" distributions in the impact parameter plane, 
27r6P(6). 



binary parton-parton collisions. Neglecting possible correlations between the partons in the 
transverse plane it is given by 

P,{h) ^ [P2{h)f I j d'h [P,{h)f- (38) 

Fig.fTSlshows that this distribution is significantly narrower than P2, i.e., the requirement of two 
hard processes results in a further reduction of effective impact parameters. 

Correlations in the transverse positions of partons can be probed by studying pp/pp events 
with two hard processes, involving two binary collisions of partons. At the Tevatron such a 
process — production of three jets and a photon — was studied by the CDF collaboration 
mm . The observed cross section is by a factor of 4 larger than the naive estimate based on 
the assumption that the partons are distributed uniformly in the transverse plane, over an area 
whose size was inferred from the proton electromagnetic form factor. The effect of correlations 
in the transverse position of partons {i.e., their localization in "spots" of significantly smaller 
size than the radius of their distribution within the nucleon) reduces this discrepancy by a factor 
of 2. This hints at the presence of significant correlations in the parton transverse positions for 
X > 0.05. A possible explanation of such correlations is the localization of the non-perturbative 
gluon fields in "constituent" quarks (and antiquarks), as suggested by the instanton vacuum 
model of chiral symmetry breaking in QCD and supported by a large body of information on 
low-energy hadron structure, see Ref. 11311 for a review. We find that the parton correlations 
implied by this model indeed give rise to a further enhancement of the cross section for two 
hard processes by a factor ~ 2, see Refs. ir73llll2il for details. 
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7.2 Black-disk limit in high-energy pA and pp collisions 



A new effect encountered in high-energy pA and pp collisions is that the interaction of lead- 
ing partons in the proton with the gluon field in the nucleus (or other proton) approaches the 
maximum strength allowed by s-channel unitarity, the BDL. This leads to certain qualitative 
modifications of the hadronic final state, which will be observable in central pp and pA colli- 
sions at LHC. In particular, this effect dramatically changes the strong interaction environment 
for new heavy particle production in central pp collisions at LHC. 

In a high-energy pA collision, consider a leading parton in the proton, with longitudinal 
momentum fraction xi ~ 10"^ and typical transverse momentum of the order of the inverse 
hadron size; see Fig. [161 Assuming a leading-twist two-body scattering process with transverse 
momentum p± in the final state, this leading parton can interact with partons in the nucleus of 
momentum fraction 

= (39) 

Xi s 

(x2 and s here refer to the effective pN collision). If s becomes sufficiently large, X2 can 
reach very small values even for sizable transverse momenta, p']_ ^ ^qcd- example, at 
LHC X2 ~ 10~^ is reached for p± ^ 2GeV. At such values of X2, the gluon density in the 
nucleus becomes large. The leading parton can be thought of as propagating through a dense 
"medium" of gluons. In this situation, the probability for the leading parton to split into several 
partons and scatter inelastically approaches unity, corresponding to the scattering from a "black" 
object. As a result the leading parton effectively (via splittings) undergoes inelastic collisions, 
losing energy and acquiring transverse momentum, until its transverse momentum is so large 
that the interaction probability becomes small, and the nucleus no longer appears "black". To 
summarize, we can say that in pA collisions the leading partons acquire transverse momenta of 
the order of the maximum transverse momentum for which their interaction with the nucleus 




Figure 16: The black-disk limit (BDL) in central pA collisions: Leading partons in the proton, 
Xi ~ 10'\ interact with a dense medium of small-a;2 gluons in the nucleus (shaded area), 
acquiring a large transverse momentum, p±. 
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Figure 17: Black-disk limit in central pA collisions at LHC: (a) The profile function for the 
scattering of a leading gluon in the proton (regarded as a constituent of a gg dipole) from 
the nucleus at zero impact parameter, r°'^(6 = 0), as a function of the transverse momentum 
squared, p\. (b) The maximum transverse momentum squared, p\ g^^' for which the interaction 



of the leading gluon is "black", F > Fcnt, as a function of the gluon's momentum fraction, xi. 
Here we assume ^/s = 14 TeV for the effective NN collisions, in order to facilitate comparison 
with the case of central pp collisions in Fig.flSl 



at their respective Xi is close to the BDL, p_l,bdl- This transverse momentum represents a new 
hard scale in high-energy hadron-hadron collisions, which appears because of the combined 
effect of the rise of the gluon density at small x and the unitarity condition.^ 

To estimate the maximum transverse momentum for interactions close to the BDL, we can 
treat the leading parton as one of the constituents of a small dipole scattering from the target. 
This "trick" allows us to apply the results of Sec.|6lto hadron-hadron scattering. In this analogy, 
the effective scale in the gluon distribution is = 4p^, corresponding to an effective dipole 
size of d ~ 3/{2p±). For simplicity, we first consider the case of central collisions of a proton 
with a large nucleus, which allows us to neglect the spatial variation of the gluon density in 
the target in the transverse direction. This amounts to approximating the transverse spatial 
distribution of gluons in the nucleus by 

GAx,p;Ql,) - ^^ii^^ ip<R^). (40) 

As a criterion for the proximity to the BDL, we require that the profile function of the dipole- 
nucleus amplitude at zero impact parameter satisfy r'^^(6 = 0) > Fcnt, see Fig. flTk . For 
an estimate, we choose Fcnt = 0.5, corresponding to a probability for inelastic interaction of 

*The kinematics of the final state produced in the interaction of the large-a;i parton with the small-X2 gluon 
field resembles the backscattering of a laser beam off a high-energy electron beam. The large-xi parton gets a 
significant transverse momentum and loses a certain fraction of its longitudinal momentum, accelerating at the 
same time a small-X2 parton. 
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Figure 18: Upper row: The transverse momentum squared, g^L, acquired by a leading gluon 
(momentum fraction xi) through interactions with the small-X2 gluon field in the other proton 
near the BDL, as a function of the impact parameter of the pp collision, h. Shown are the 
estimates for LHC (left panel) and Tevatron energies (right panel). Lower row: Average values 
of p^ in pp collisions with a single hard process (impact parameter distribution P2) and two 
hard processes (distribution P4), cf. Fig.[T5l For leading quarks, the values of {p\_ b^l) are about 
half of those for gluons shown here. 



0.75, reasonably close to unity. We then determine the maximum px_ for which the criterion 
is satisfied. Fig. [l71 shows the result for p\ for a leading gluon, as a function of the gluon 
momentum fraction, xi \ for leading quarks, the result for p\ g^L is approximately 0.5 times the 
value for gluons. The numerical estimates show that leading partons indeed receive substantial 
transverse momenta when traversing the small-X2 gluon medium of the nucleus. We emphasize 
that our estimate of p±,bdl applies equally well to the interaction of leading partons in the 
central region of A A collisions. 

Turning now to pp collisions, we have to take into account the transverse spatial structure 
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of the colliding hadrons. A crucial point is that high-energy interactions do not significantly 
change the transverse position of the leading partons, so that their interaction with the small- 
X2 gluons is primarily determined by the gluon density at this transverse position. Because 
the leading partons in the "projectile" proton are concentrated in a small transverse area, and 
the small-a;2 gluon density in the "target" proton decreases with transverse distance from the 
center, it is clear that the maximum transverse momentum for interactions close to the BDL, 
pi BDL' decreases with the impact parameter of the pp collision, b. Fig. [TS] (upper row) shows 
the dependence of p^ obtained with the parametrization of the transverse spatial 

distribution of gluons based on analysis of the HERA exclusive data, see Section|43||73 |. One 
sees that p_l,bdl ~ several GeV in central collisions at LHC. Substantially smaller values are 
obtained at the Tevatron energy. 

To determine the typical transverse momenta of leading partons in events with new particle 
(or hard dijet) production, we need to average the results for over pp impact parameters, 

with the distribution implied by the hard production process, P2{b), Eq. (l37t . or, in the case of 
four jet production, with P^i^b), Eq. (l38b . The resulting average values of p^ are shown in 
Fig. HE] (lower row). We find that the suppression of large impact parameters implied by the 
hard process, described in Sec. 17. 11 is sufficient to keep p_l,bdl above 1 GeV in more than 99% 
of events at LHC. 

To summarize, our estimates show that in generic central pA and central (triggered) pp col- 
lisions at LHC the leading partons acquire substantial transverse momenta due to interactions 
near the BDL. A much weaker effect is found at the Tevatron energy. The origin of this differ- 
ence is the increase in the gluon density due to the decrease of X2 between Tevatron and LHC 
energies, cf. Eq. (l39l) . 

7.3 Final state properties in central pp collisions at LHC 

The approach to the BDL in the interaction of leading partons implies certain qualitative chan- 
ges in the hadronic final state in central pp and pA collisions. In particular, these effects will 
profoundly influence the strong interaction environment for the production of new heavy parti- 
cles (Higgs boson, etc.) at LHC. 

The main effect of the BDL is that the leading partons in the projectile acquire substan- 
tial transverse momenta, of the order when propagating through the dense medium of 
small-a;2 gluons in the target. As a result, the projectile becomes "shattered": The leading 
partons lose coherence and fragment independently over a wide range of rapidities close to the 
maximal rapidity, corresponding to hadron momentum fractions z ~ xiyu/p^ bdl (/U is a typical 
hadronic mass scale). The differential multiplicity of leading hadrons, integrated overpx, is ap- 
proximately given by the convolution of the nucleon parton density, /„, with the corresponding 
parton fragmentation function, Dh/a, at the scale Qgff = 4p5_ bdl' 



where is total number of inelastic events II113LI1061I114II . This corresponds to a very strong 
suppression of forward hadron production as compared to generic inelastic pp events. The sup- 
pression is particularly pronounced for nucleons; one expects that for ^ > 0.1 the differential 




(41) 
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multiplicity of pions should exceed that of nucleons. At the same time the transverse spectrum 
of the leading hadrons will be much broader, extending up to p±^bdl ^ 1 GeV. Finally, the 
independent fragmentation mechanism implies that there will be no correlations between the 
transverse momenta of leading hadrons (some correlations will remain, however, because two 
partons produced in collisions of large-xi and small-a;2 partons may end up at similar rapidi- 
ties). 

In central pp collisions at LHC, where leading particles are suppressed in both the forward 
and backward direction, one expects a large fraction of events with no particles with z > 0.02 — 
0.05 in both fragmentation regions. This amounts to the appearance of long-range rapidity 
correlations. Such events should show a large energy release at rapidities y = A — 6. However, 
similar to the case of diffractive processes in ep scattering (cf. the discussion in SectionlU), one 
should expect that there is a ~ 10% probability for dijets to be produced in pp collisions at large 
impact parameters without additional interactions between the constituents of the nucleons. 

Another important effect of the BDL is a significantly increased energy loss of the leading 
partons, due to the larger probability of inelastic collisions, and the wider distribution of the 
propagating parton (dipole) over transverse momenta. In particular, a 2% energy loss would 
explain the suppression of forward pion production in deuteron-nucleus collisions at RHIC 
at pj_ ~ 4 GeV f99\. Studies of this effect would be possible both at RHIC and LHC, in 
particular if the forward capabilities of the current detectors were upgraded, as discussed in 
several proposals presently under consideration. Note that energy loss is neglected in Eq. (|4TT) . 
This corresponds to the usual assumption of models in which parton propagation is treated 
as multiple elastic rescattering of the parton's accompanying gluon field from the medium. A 
consistent treatment of energy loss and transverse momentum broadening near the BDL remains 
a challenge for theory. We note that the pattern of energy loss in our approach is qualitatively 
different from models in which the leading partons scatter from a classical gluon field (in that 
case energy loss is negligible) HllSII . 

The approach to the BDL has consequences also for hadron production in the central rapidity 
region. The multiple scattering of large-xi projectile partons from the small-X2 gluons in 
the target shifts a large number of these gluons to larger rapidities, leaving numerous "holes" 
in the target wave function. Furthermore, multiple interactions of partons with moderately 
small xi ~ X2 also occur with large probability. (Unitarity effects should be important for 
these interactions as well, but have not been studied so far.) Both effects lead to the creation 
of a substantial amount of color charge, which should result in an increase of soft particle 
multiplicities over a broad range of rapidities as compared to the situation far from the BDL. 
This increase should in fact be observable already at Tevatron energies, in central events selected 
by a trigger on two-jet or production. An increase of the multiplicity at rapidities \y\ < 1.0 
in such events compared to minimum bias events was indeed reported in Ref. II116II . It would 
be extremely interesting to extend such studies to higher rapidities. 

Our findings imply that new heavy particles at LHC will be produced in a much more "vio- 
lent" strong-interaction environment than one would expect from the extrapolation of the prop- 
erties of minimum bias events at the Tevatron. Even the extrapolation of properties of hard 
dijet events should not be smooth, as the transverse momenta acquired by leading partons are 
estimated to be substantially larger at LHC than at Tevatron, see Fig.fTSl 
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7.4 Black-disk limit in elastic pp scattering 



The assumption of the BDL in the interaction of leading partons, combined with the complex 
structure of the proton wave function in QCD, allows us to estimate the profile function of the pp 
elastic amplitude at small impact parameters. This simple estimate nicely explains the observed 
"blackness" of the phenomenological pp profile function at 6 = at the Tevatron energy, and 
allows us to extrapolate the profile function at small b to higher energies. It also raises the 
question whether the observed blackness could be explained on the basis of soft interactions. 

In pp collisions at small impact parameters, leading quarks on average receive substantial 
transverse momenta when passing through the small-x gluon field of the other proton; see 
Fig. [HI When a single leading quark gets a transverse momentum, p±, the probability for 
the nucleon to remain intact is approximately given by the square of the nucleon form factor, 
F^{p\), which is < 0.1 for p^_ > 1 GeV, i.e., very small. One may thus conclude that the 
probability of the survival averaged over should be less than 1/2 (on average, half of the time 
the quark should receive a transverse momentum larger than the average transverse momentum, 
p± > I GeV). Because there are six leading quarks (plus a number of leading gluons), the 
probability for the two protons to stay intact, |1 — rPP(s,6)p, cf. Eq. (l29b . should go as a 
high power of the survival probability for the case of single parton removal, and thus be very 
small. This crude estimate shows that already at Tevatron energies |1 — r^P(s, 6 = 0)p should 
be close to zero owing to hard interactions. The conclusion that the small impact parameter 
hadron-hadron interactions should become black at high energies follows principally from the 
composite structure of the hadrons and does not depend on any details. In particular, if taming 
effects stopped the growth of the dipole-nucleon interaction at a fraction of the BDL, our result 
would not change. Our conclusion that TP^{y/s > 2TeV, 6 = 0) ~ 1 agrees well with the 
current analysis of the Tevatron data, see e.g. Ref. ifTTOI . 

One can estimate the maximum impact parameter, 6^7, up to which hard interactions cause 
the pp interaction to be "black". The probability for a leading parton with xi ~ 10~^ to expe- 
rience a hard inelastic interaction increases with the collision energy at least as fast as dictated 
by the increase of the gluon density in the other proton at X2 = 'ip']_/{sxi), cf. Eq. (l39l) be- 
low. Because X2G{x2,Q'^) oc Xa"^'' with Uh > 0.2, the probability should grow as s"''.^ The 
dipole parametrization of the transverse spatial distribution of gluons, Eq. (l22t . suggests that 
the gluon density decreases with the distance from the center of the nucleon approximately as 
~ exp [—mg(x2)p]. If one neglects the transverse spread of the large-xi partons as compared 
to that of the small-X2 gluons one arrives at an estimate of the energy dependence of 6j? as due 
to hard interactions II102I . 

i,. « !!iMfM, (42) 

mg{x2) 

where ^/sr = 2 TeV is the Tevatron energy. In principle, may decrease at very large vir- 
tualities, which would become important at extremely high energies. However, this effect is 
likely to be compensated by the increased number of constituents in the nucleon wave function 
affected by the BDL. The above estimate is consistent with the popular Pomeron model param- 
eterization of the pp elastic amplitude II117II . In this model Tpp(0) ^ 1 for s = st and s-channel 

^The HERA data on dipole-nucleon scattering suggest that the taming of the gluon density starts only when 
the probability of inelastic interactions becomes large, > 1/2. However, for such probabilities of single parton 
interactions, multiparton interactions ensure that the overall interaction is practically black. 
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Figure 19: Energy dependence of the maximum impact parameter for "black" pp interactions, 
hp- Solid line: Estimate based on the pomeron model fit of pp elastic scattering of [.117.1 . Dashed 
line: Estimate based on hard interactions, Eq. (H^ . 

unitarity is violated dXh < hp for s > st- The dependence of 6^? on the energy in the Pomeron 
model is similar to our estimate (l42t . see Fig. [191 This shows that the two-scale picture of the 
transverse structure of the nucleon — and the ensuing picture of hard and soft interactions — 
are self-consistent. 

Overall, the analysis of Ref. ["1021 summarized in this subsection suggests that at high en- 
ergies hadron-hadron interactions should be "black" in a range of impact parameters growing 
approximately like In s, with a coefficient growing only very slowly with energy due to de- 
crease of mg{x2) with energy. This corresponds to the Froissart regime, with interactions in the 
black region dominated by semi-hard interactions, and interactions at large impact parameters 
dominated by the single Pomeron exchange. At the same time, this analysis indicates that the 
concept of summing multi-Pomeron exchanges, which should give the dominant contribution 
at small b, breaks down as the soft physics is gradually squeezed out to large b. 

It is worth emphasizing here that, in principle, the BDL can emerge in hadron-hadron in- 
teractions already at the level of the soft interactions 111181 fTT 9l. However, it is hardly possible 
to reconcile it with the pre-QCD Feynman parton model description of high energy processes, 
if one would require it to be valid both in the rest frame of the target and in the center-of- 
mass frame. Really, within the parton model one cannot generate complete absorption of the 
projectile in central collisions in the target rest frame, where the target consists of few partons 
lfT20J . This puzzle is resolved in QCD, where radiation leads to the "blackening" of the hard 
interactions at central impact parameters. At extremely high energies, as a result of this effect, 
all memory of the colliding hadrons is lost. Hence the universal behavior of total cross sections, 
(Ttot cx; ln^(s/so), with universal coefficient for all hadrons (nuclei) II102I . ^ 

^The increase of the interaction with energy, and the related increase of essential impact parameters, show 
that the theoretical description of high-energy hadronic collisions should be closer to classical mechanics than to 
quantum mechanics (V. Gribov, private communication to Yu. Dokshitzer). An example is the cross section for the 
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7.5 Ion-induced quark-gluon implosion 



The small-x phenomena outlined above — the approach to the BDL, and large leading-twist 
gluon shadowing — play an important role also in the heavy ion collisions at LHC energies. 
Here we consider just one example, the so-called ion-induced quark-gluon implosion in the 
nucleus fragmentation region. For a review of other effects in the framework of the color glass 
condensate model, see Ref. II115I . 

In generic central AA collisions at collider energies, in analogy to the central pA collisions 
discussed in Section fT^ all the leading partons of the individual nucleons are stripped off "soft" 
partons and form a collection of quarks and gluons with large pj_. In the rest frame of the frag- 
menting nucleus, the incoming nucleus has a "pancake" shape with longitudinal length ~ 1 fm 
for soft partons, and _Rjv("^Af/PAr)(xy/x) <^ for hard partons, where xy ~ 0.2 is the av- 
erage X value for the valence quarks. That is, the nucleons in the nucleus at rest at different 
locations along the collision axis are hit by the hard partons in the incoming nucleus one after 
another. In the BDL, no spectators are left. The hit partons are produced with practically the 
same x that they had in the nucleus (because the fractional energy loss is small), transverse mo- 
menta ~ p±,BDL, and virtualities < p^ g^L- The partons move in the direction of the projectile 
nucleus. Because they are emitted at finite angles, their longitudinal velocity is smaller than 
the speed of light, and they are left behind the projectile wave. However, because the emission 
angles are small, a shock wave is formed, compressing the produced system in the nucleus rest 
frame. In the frame co-moving with the shock wave, valence quarks and gluons are streaming 
in the opposite directions. The resulting pattern of fragmentation of the colliding nuclei leads 
to an "implosion" of the quark and gluon constituents of the nuclei. The non-equilibrium state 
produced at the initial stage in the nucleus fragmentation region is estimated to have densities 
oc pi BDL' which is > 50GeV/fm^ at LHC, and probably > lOGeV/fm^ at RHIC. It seems 
likely that the partons would rescatter strongly at the second stage, although much more de- 
tailed modeling is required to find out whether the system would reach thermal equilibrium. 
Such large-angle rescattering of partons would lead to production of partons at higher rapidi- 
ties, and re-population of the cool region. In particular, two gluons from the pancake could 
have the right energies to produce near-threshold cc pairs and Xc-mesons with small transverse 
momenta and xf{cc) ~ 2xg ~ O.L 

7.6 Cosmic ray physics near the GZK cutoff 

An extensive program of studies of cosmic rays at energies close to the Greisen-Zatsepin- 
Kuzmin (GZK) cutoff [TO"], Eqzk — 6 x 10^° GeV, is under way, using several cosmic ray 
detectors. These experiments detect cosmic rays indirectly, via the air showers induced when 
they enter the atmosphere. The properties of the primary particle need to be inferred from 
those of the observed shower. For this, a good understanding of the physics of high-energy 
interactions in the atmosphere is mandatory. The observed characteristics of the shower are 
predominantly sensitive to leading hadron production (xp > 10^^), which, according to our 

scattering of a high-energy particle from a potential rapidly decreasing with impact parameter The essential impact 
parameters — and therefore the cross section — are infinite within classical mechanics, but finite in quantum 
mechanics, while they increase with energy in QCD. 
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discussion above, at these energies probes small-x dynamics down to x ~ 10^ , deep inside 
the regime affected by the approach to the BDL. First studies of these effects were performed in 
Ref . [ 1 2 1 1 . It was found that the steeper x ^^-distribution of leading hadrons as compared to low- 
energy collisions, caused by the strong increase of the gluon densities at small x (see Section 
YlAk leads to a reduction of the position of the shower maximum, X^ax- Account of this effect 
in the models currently used for the interpretation of the data may shift fits of the composition 
of the cosmic ray spectrum near the GZK cutoff towards lighter elements. Also, it appears that 
the present data on X^^^^E) exclude the possibility that the prediction of a rapid growth of the 
critical x-value where the BDL becomes effective (~ which is compatible with RHIC 

and HERA data, would persist up to the GZK cutoff energy. 

8 Hard diffraction at hadron colliders 
8.1 Diffractive proton dissociation into three jets 

LHC will offer an opportunity to study a variety of hard diffractive processes in pp and pA 
scattering. One interesting aspect of such processes is that they allow to probe rare small- 
size configurations in the nucleon wave function.A proton in such a configuration can scatter 
elastically off the target and fragment into three jets, corresponding to the process 

p + p{A) jetl +jet2 +jet3 (43) 

The cross section for the diffractive process can be evaluated based on the kind of QCD 
factorization theorem derived in Ref. [ 19J . The cross section is proportional to the square of the 
gluon density in the nucleon at a; ~ M^(3 jets)/s, and virtuality ~ (1 — 2)p5_ II122I . The 
distribution over the fractions of the proton longitudinal momentum carried by the jets is propor- 
tional to the square of the light-cone wave function of the \qqq) configuration. The numerical 
estimates suggest II123I that the process could be observed at the LHC energies provided one 
would be able to measure jets withp^ ~ 10 GeV at very high rapidities, ?/jet(p_L = lOGeV) ~ 6, 
and with a large background from leading-twist hard diffraction. The latter will be suppressed 
in pA collisions, because the coherent 3-jet process has a much stronger A-dependence than 
the background due to soft and hard diffraction induced by strong interactions. The main back- 
ground will be due to hard electromagnetic interactions of the proton with the Coulomb field 
of the nucleus. We note that the discussed mechanism of hard diffraction requires that the in- 
teraction of the spatially small three-quark color singlet configuration with the proton be far 
from the BDL at small impact parameters. Otherwise production at small impact parameters 
would be suppressed, leading to a dip in the t-dependence of the differential cross section for 
the production of three jets with moderate p±. 

The detection of the three-jet final state produced by diffractive scattering of a qqq config- 
uration from a proton should be easier than that resulting from e+e^ annihilation into qqg, as 
in the former case all color charges are in the triplet representation, leading to less radiation be- 
tween the jets. Finally, it would also be possible to study the process pp ^ pn + two jets, which 
is similar to pion dissociation into two jets. Experimentally, this would require the measurement 
of jets at rapidities y ~ 4, together with the detection of a leading neutron by a zero-degree 
calorimeter, as is present in several of the LHC detectors (ALICE, ATLAS, CMS). 
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8.2 Exclusive diffractive Higgs production 



Hard diffractive processes are also being considered in connection with the production of new 
heavy particles in pp collisions at LHC. In particular, the exclusive diffractive production of 
Higgs bosons, 

p + p p + (rapidity gap) + H + (rapidity gap) + p, (44) 

is regarded as a promising candidate for the Higgs search; see Ref. 11241 and references therein. 
From the point of view of strong interactions, this process involves a delicate interplay between 
"hard" and "soft" interactions, which can be described within our two-scale picture of the trans- 
verse structure of the nucleon II112II . The Higgs boson is produced in a hard partonic process, 
involving the exchange of two hard gluons between the nucleons. The impact parameter dis- 
tribution of the cross section for this process is described by the square of the convolution of 
the transverse spatial distributions of gluons in the in and out states, P4(6), defined in Eq. (l38t . 
where the scale is of the order of the gluon transverse momentum squared, ~ Mfj/A. In ad- 
dition, the soft interactions between the the spectator systems have to conspire in such a way 
as not to fill the rapidity gaps left open by the hard process. The probability for this to happen 
is approximately given by one minus the probability of an inelastic pp interaction at a given 
impact parameter, or |1 — r^^(s, 6) p. The product of the two probabilities, which determines 
the 6-distribution for the total process, is shown in Fig.l20h. At small h the probability for no 
inelastic interaction is very small |1 — F^^p ^ 0, leading to a strong suppression of small h in 
the overall distribution. 

The so-called rapidity gap survival probability, which measures the "price" to be paid for 
leaving the protons intact, is given by the integral HI 1211 

5' = J d%\l-TPP{s,b)\^ Piib). (45) 

Fig. EDb shows our result for this quantity, with s ranging between Tevatron and LHC ener- 
gies, for various values of the dipole mass in the two-gluon form factor of the nucleon, m^, 
Eq. (l22b . The survival probability decreases with s because the size of the "black" region at 
small impact parameters (in which inelastic interactions happen with high probability) grows 
with the collision energy. Note that the effective x values in the gluon distribution decrease with 
the energy (for fixed mass of the produced Higgs boson), resulting in smaller effective values 
of m^. This makes the actual drop of the survival probability with energy slower than appears 
from the fixed-m^ curves of Fig. l20b . Our estimates of S*^ are in reasonable agreement with 
those obtained by Khoze et al. II125I in a multi-Pomeron model, as well as with those reported 
by Maor et al. [il26j . In view of the different theoretical input to these approaches this is very 
encouraging. 

Our results for the rapidity gap survival probability apply equally well to the production 
of two hard dijets instead of a Higgs boson. For this process, one expects much larger cross 
section, and it would be possible to investigate experimentally the interplay of hard physics and 
absorptive effects, which leads to a rich, distinctive structure of the cross section as a function 
of the transverse momenta of the two protons, Ai^ and 111271 . This structure should also 
rather strongly depend on the rapidities of the jets, due to the x-dependence of the transverse 
spatial distribution of gluons, see Sec. 14.41 (L. Frankfurt et al., in preparation). 
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Figure 20: (a) The impact parameter (b-) distribution of the cross section for diffractive Higgs 
production at LHC (-\/s = 14TeV). Dashed line: 6-distribution of the hard process, ^4(6), 
Eq. (|38l), cf. Fig.ini Solid line: 6-distribution of the total process, |1 - V'p'p{s, h)\'^Pi{h). Here, 
m?g = 1 GeV^. (b) The rapidity gap survival probability, 5^, Eq. (l45l 1 1 12|. Shown is the result 
as a function of s, for various values of the mass parameter in the two-gluon form factor, m^. 
The Tevatron and LHC energies are marked by arrows. 



8.3 Inclusive hard diffractive processes 



Inclusive hard diffractive processes, such as 



p + p p + (rapidity gap) + 2 jets + X, 
p + p — > p + (rapidity gap) + 2 jets + X 



(rapidity gap) 



(46) 



P, 



offer a possibility to probe the "periphery" of the proton with hard scattering processes. The 
cross section for these processes is again suppressed compared to the naive estimate based on 
the diffractive parton densities of the proton measured in ep scattering at HERA. As in the case 
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of exclusive diffractive Higgs production, the cause of this is the very small probability for the 
nucleons not to interact inelastically at small impact parameters. The suppression factors can be 
estimated by generalizing the approach to the description of hard and soft interactions outlined 
in Sec. 18.21 Simple estimates along the lines of Eq. (l45t naturally reproduce the suppression 
factors of the order 0.1 — 0.2 observed at Tevatron. However, the results in this case are more 
sensitive to the details of the impact parameter dependence of the hard scattering process and 
the soft spectator interactions. 

Measurements of inclusive hard diffractive processes at LHC would allow one to perform 
many interesting tests of the diffractive reaction mechanism. In particular, one could (a) inves- 
tigate how the overall increase of the nucleon size with energy leads to a suppression of hard 
diffraction, (b) check how the rate of suppression depends on the x-value of the parton involved 
in the hard process, (c) look for the breakdown of Regge factorization, that is, the change of the 
diffractive parton distributions with xp. 



9 Summary and Outlook 
9.1 From HERA to LHC 

The HERA experiments and the theoretical investigations they stimulated have greatly advanced 
our knowledge of small-x dynamics. The key result of these studies are: (a) The rapid increase 
with energy of the cross section for the scattering of small qq wave packets from the nucleon. 
HERA energies are not sufficient to reach the BDL in the dipole-nucleon interaction in average 
configurations. The interaction of gluon dipoles in diffractive scattering appears to be close to 
the unitarity limit for ~ 4 GeV^, but this can hardly be verified directly because of the lack 
of a trigger for such configurations, (b) The establishment of a three-dimensional picture of 
the partonic structure of the nucleon. The leading partons are concentrated in a much smaller 
transverse area than the area associated with the nucleon in soft hadronic processes at high 
energies. 

We have demonstrated that these elements of small-a; dynamics are of utmost importance 
for building a realistic description of pp/pA collisions at LHC. The BDL will be commonplace 
in central pp/pA collisions at LHC, affecting average configurations in the colliding protons 
(nuclei), with numerous consequences for the hadronic final state. In particular, these phenom- 
ena qualitatively change the strong interaction environment for new particle production. 

We have identified several directions for future theoretical research, necessary for describing 
the expected new phenomena at LHC. These include the resummation approaches to QCD 
radiation (combining logarithms of and the account for energy loss in the interaction 
of leading partons with the small-x gluon medium, and the development of realistic models of 
hadron production in central pp/pA collisions with interactions close to the BDL. 

An overarching goal of future theoretical research on the structure of the nucleon and small- 
X dynamics should be to bring together the approaches starting from "soft" (hadronic) and 
"hard" (partonic) physics, as envisioned in Gribov's space-time picture of high-energy inter- 
actions. We have pointed out several instances in which "soft" and "hard" dynamics match 
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smoothly or exhibit a delicate interplay, e.g. pp elastic scattering at small impact parameters, or 
diffractive processes. 

A natural question is what are the most promising directions for future experimental studies 
of small-x dynamics. This question needs to be discussed with regard to the general long- 
term perspectives in high-energy physics. We assume here that the decision by the DESY 
management to stop HERA operations in 2007 will be enacted. This would clearly be a great 
loss, as many insights could be obtained from further measurements at HERA in both ep and 
eA mode, see e.g. the proposals put forward for the HERA III run. We shall thus focus on the 
possibilities offered by LHC, with comments on the possible future program at Tevatron, as well 
as on the electron-ion collider envisaged in the U.S. government's long-range plan 1 128 1. The 
small-x investigations at LHC described in the following subsections are meant to complement 
the studies of central inelastic pp/pA collisions (Sec. 17 .41) and diffractive phenomena in pp/pA 
scattering, which are the main topic of this review. 

9.2 Measurement of parton densities in pp and pA collisions at LHC 

Measurements at LHC could greatly expand the x-range in which the parton densities are 
known. This would require measurements of hard processes such as 



in the region where one of the colliding partons carries small momentum fraction. The cross 
sections of all these processes remain large down to the very edge of the LHC kinematics, 
corresponding to x ~ 3 x 10~^ for Drell-Yan pair production with M^+^- = 5 GeV [1291 
I123II . The main limitations come from the need to identify relatively low-px jets, and from the 
detector acceptance. The smaller the x one wants to probe, the more forward one must look, as 
the momentum fractions of colliding partons are related to the rapidities of the produced jets as 



The presently planned configuration of the CMS detector would allow for the measurement of 
dijet production down to x ~ 3 x 10^^ at p± = 10 GeV. This would push parton distribution 
measurements deep into the region where unitarity effects play an important role in the dynam- 
ics of hard processes, and where evolution effect in both ln(l/x) and InQ^ need to be taken 
into account. When the BDL is reached, the Af^ -dependence of the cross section, da / dxidx2, 
is predicted to be much slower than oc 1/IVP as in the leading twist approximation, similarly to 
the case of the inclusive deep inelastic scattering II130II . 

If several of the reactions (l47t were measured, it would allow for independent tests of the 
QCD factorization, which may be violated at intermediate virtualities owing to the strong inter- 
action of the propagating system with the small-x gluon medium, see Sec.|71 The latter will be 
strongly enhanced in the region of xi close to 1. These effects, which are of great interest in 
themselves, can be probed by comparing the production cross sections for fixed, large xi and 
various values of X2, including relatively large ones where the parton densities are known. 



jetl + jet2 + X dijet production 

jet + 7 + X, 7 + 7 + X photon production 

Q + Q + X heavy quark production 

+ l~ + X Drell-Yan pair production 

W^{Z) + X weak boson production, 



(47) 




(48) 
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9.3 Small-x phenomena in ultraperipheral collisions at LHC 



It has been long known that nuclei in high-energy collisions generate a large flux of equivalent 
photons, which are spread in the transverse plane over distances substantially larger than twice 
the nuclear radius — the maximal distance at which strong interactions are possible. Scattering 
processes induced by these photons are referred to as ultraperipheral collisions. They have a 
distinctive signature, which allows them to be separated from the more frequent events caused 
by strong interactions. Experimentally, one selects events in which one of the nuclei remains 
intact, or emits one or a few neutrons by way of dipole excitation. Such events are extremely 
rare in scattering at impact parameters smaller than twice the nuclear radius. 

The experiments at HERA have shown that photon-induced processes provide a well- 
understood probe of the gluon density in the proton. At LHC, such processes could be studied 
up to invariant 7p energies (i.e., 'jA energies per nucleon) exceeding the maximal HERA energy 
by a factor of 10. This would allow one to use dijet (charm, etc.) production to measure the 
gluon density in the proton/nucleus down to s ^ 3 x 10~^ for p± ~ few GeV FlBUIlBll . as 
well as the diffractive gluon density. Among other things, measurements of diffractive chan- 
nels would allow one to perform critical tests of the HERA observation of a large probability of 
gluon-induced diffraction (see Sec.|51), and the prediction of its further enhancement for nuclear 
targets, see Ref. HI 091 and references therein. Another important measurement would be the t- 
dependence of gluon-induced diffraction and its change with energy, using CMS-TOTEM in 
pA mode. We remind the reader that the lack of direct information on the t-dependence of 
diffraction leads to a large uncertainty in the predictions for leading-twist nuclear shadowing 
(see Sec. 15.31) . 

Ultraperipheral collisions would also allow one to study the coherent production of heavy 
quarkonia, 7^4 J/tp (T) + Aatx < 10"^, and to investigate the propagation of small dipoles 
through the nuclear medium at high energies. The QCD factorization theorem predicts that the 
A-dependence of the amplitude for this process should change between the color transparency 
regime (observed at FNAL [1771 ). where it is oc A, and the perturbative color opacity regime, 
where it is proportional to the leading-twist shadowed gluon density. It would be possible also 
to use coherent diffraction from nuclei to study the approach to the BDL in 'jA — > X + A, by 
comparing the measured cross section to the BDL prediction, cf. Sec. 16.41 The most promising 
channels are J/'ip and dijet photoproduction; see Ref. II133I for a review and discussion. In 
AA collisions, it is difficult to separate processes induced by the photons generated by the 
left- and right-moving nucleus. Away from zero rapidity, a low-energy contribution tends to 
dominate, limiting the range of x which could be explored for production of a state with mass 
M to X > MA/{2pa), where pa is the momentum of the colliding nuclei. However, it seems 
that selection of events in which the heavy nucleus undergoes a dipole excitation enhances the 
contribution of hard photons [|134I . allowing one to extend the x-range of the measurements 
(by a factor of up to 10 in the case of J/^ production). The challenge is to trigger both on 
events with and without break-up. In pA mode, the dominant process will be the production of 
heavy quarkonia. Such measurements would extend the VT-range of the HERA measurements 
by a factor of three, and make it possible to measure directly the t-dependence of the cross 
section in a very broad range of rapidities, using the proposed 420m proton tagger II135I . which 
is critical for a more accurate determination of the x-dependence of the nucleon's transverse 
structure, see Sec. 14.41 Note also that in pp scattering it is possible to detect protons at very 
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small momentum transfers, where Coulomb exchange dominates II136I . This would allow one 
to measure exclusive photoproduction of heavy quarkonia in pp scattering with good statistics 

9.4 Small-x physics at RHIC and an electron-ion collider 

The LHC measurements described in Sections W?^ and 1^31 will probe small-x dynamics at least 
down to X ~ 10^^. However, most of these measurements are restricted to scales Q >h GeV, 
and it would be difficult to connect them with the physics at smaller scales (virtualities), Q ~ 
2 — 3 GeV, relevant for the overall structure of central pp/pA collisions at LHC (see Sec. |7l). 
The gap could be filled, to some extent, by experiments at RHIC and the proposed electron-ion 
collider lfT28lfT39ll . 

Extension of the forward acceptance of the current RHIC detectors would make it possible to 
measure Drell-Yan pair production at x ~ 10"'^ in pp and pA/dA scattering. This would allow 
one to test the predictions for leading-twist nuclear shadowing and look for deviations from 
the leading-twist prediction in the pi_ distributions of the dileptons. Qualitatively, one expects 
a suppression of the low transverse momentum part of the distribution up to p±^ ~ bdl- 
As mentioned above, such measurements would also allow one to probe the role of final state 
interactions by varying the x-values of the leading partons in the proton. If absorption effects 
were significant, one would have to introduce a cut on Xp < 0.3 to suppress these effects, which 
would reduce somewhat the x^-range where the parton densities can be probed. 

The eRHIC design for a future electron-ion collider envisages an ep/eA collider with 
-y/i < 100 GeV, with significantly higher luminosity than HERA, and the ability to continu- 
ously vary the beam energies over a wide range 11281 11391 . With such a facility one could 
systematically study a variety of color transparency phenomena and use them to disentangle 
the quark-gluon structure of hadrons and nuclei; one could also measure longitudinal cross 
sections, which provide stringent tests of the range of the validity of the leading-twist approx- 
imation at small x (see Sections |3l |4| and |5l). In eA collisions, one could study the transition 
of the nonperturbative shadowing at low Q"^ to the regime of leading-twist shadowing at high 
Q^, and explore whether there exists an "intermediate" regime characterized by weak coupling 
but large parton densities. The ability to perform such measurements with a range of nuclear 
beams would allow one to study these effects as a function of the nuclear thickness, reaching 
values 1 .5 times larger than the average thickness of the heavy nuclei. No other planned fa- 
cilities would be able to cover this important kinematic region. Finally, eRHIC would make it 
possible to measure the t-dependence of a variety of hard exclusive processes in a wide range 
of X, 0.1 > a; > 0.003. This would probe the transverse structure of the proton directly in the 
x-range relevant for understanding nucleon fragmentation in central pp/pA collisions at LHC. 
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